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ABSTRACT 


Ocean  currents  and  surface  waves  may  induce  serious 
errors  in  oceanographic  measurements  obtained  from  moored 
and  drifting  buoy  systems.  A general,  computationally 
efficient  solution  for  the  dynamics  of  moored  buoy  systems, 
free  drifting  buoys,  and  drogued  buoy  systems  in  three- 
dimensional  space  is  described  in  this  report.  Time-domain 
computer  simulations  of  four  specific  configurations  in 
various  environments  are  presented.  The  mathematical  model 
of  one  of  these  configurations,  a subsurface  mooring,  was 
evaluated  and  improved  using  full  scale  ocean  test  data. 

The  model  of  the  surface  mooring  configuration  studied  will 
soon  also  be  evaluated  using  recent  test  data. 
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NOMENCLATURE 


Acceleration  vector. 


AN'AT 


Area  used  in  drag  calculations. 
Acceleration  of  point  c in  the  E frame. 

Normal  and  Tangential  areas  for  a 
cylindrical  body. 

Relative  acceleration  vector. 


Point  of  attachment  between  the  surface 
buoy  and  the  mooring  line. 

Added  mass  force  vector. 


c,x‘ ,y ' ,z' 


CN'CT 


CDN'CDT 


CDIN,CDIT,  CDIA 


Body  coordinate  system,  Frame  B. 

Mid  point  of  the  mooring  line  element. 
Array  of  constant  nodal  forces. 

Drag  coefficient. 

Added  mass  coefficients,  normal  and 
tangential  directions,  surface  buoy. 

Normal  and  tangential  drag  coefficients 
for  a cylindrical  body. 

Pressure  drag  coefficient  for  the  end 
plate  of  a spar  buoy. 

Appropriate  drag  constants  corresponding 
to  DN,  DT,  and  DA. 

Viscous  drag  force  on  a body  inserted  in 
a mooring  line. 

Viscous  drag  force  vector. 

Moment  vector  due  to  viscous  drag  force. 

Normal  and  tangential  drag  forces  per  unit 
stretched  length  of  the  mooring  line. 
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dQ 


dM 


{FA} 

FF 

{FG} 

FH 


g 

H 


h 


o 


hi 


IF 


IFb 

Jxx' Iyy 


AAA 


AAA 


Viscous  drag  forces  acting  on  the 
node  of  the  mooring  line. 

Mass  of  water  displaced  by  a differential 
disk  of  the  surface  buoy,  . S dz'. 

Mass  of  a differential  disk  of  the  surface 
buoy . 

Total  force  vector  acting  on  the  body. 

Array  of  additional  forces  and  moments. 

Array  of  wave  exciting  forces  and  moments. 

Froude  Krylov  exciting  force  vector. 

Array  of  hydrostatic  and  gravitational 
forces  and  moments. 

Hydrostatic  pressure  force  acting  on  the 
mooring  line,  per  unit  stretched  length. 

Acceleration  due  to  gravity. 

Length  of  the  surface  buoy. 

Instantaneous  draft  of  the  surface  buoy, 
measured  up  to  the  mean  free  surface. 

Spar  buoy  drum  height. 

Spar  buoy  mast  height. 

hQ  plus  the  wave  elevation  component. 

Inertia  force  per  unit  stretched  length 
of  the  mooring  line. 

Inertia  force  vector  of  the  inserted  body. 

Moments  of  inertia  of  the  spar  buoy  in 
roll  and  pitch  about  c. 

Unit  vectors  of  o,x,y,z. 

Unit  vectors  of 


Wave  number  vector. 


Wave  number  component,  ;K.  cos  r . 

Wave  number  component,  jK  sin  f. 

Stiffness  coefficient  of  the  n*-^  segment. 

Integration  limits  as  defined. 

Total  mass  of  the  spar  buoy. 

Matrix  of  tensors  of  inertia  and  added 
inertia. 

Added  mass  of  the  inserted  body,  shape 
other  than  cylindrical. 

Mass  of  the  inserted  body. 

t*  H 

Mass  of  water  displaced  by  the  n'  mass. 

Added  mass  components  (normal  and 
tangential)  due  to  cylindrical  bodies. 

Number  of  nodes  or  segments  of  a mooring 
line. 

Net  bouyancy  of  the  inserted  body. 
Arbitrary  line  segment,  or  lumped  moss. 
Fixed  coordinate  system,  frame  E. 
Hydrostatic  pressure  of  the  fluid. 


Location  of  the  centroid  of  the  differential 
disk  of  a surface  buoy. 

= J" 2 dBeKZw (Zp) 1 

Displacement  vector. 

Displacement  vector  from  point  o to  point  c. 
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*'o 


.S' 


s 
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Vc 


TM 

TF 
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V 

UR 
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02 


of  t ho  spar  buoy  or  the  reduced 
radius  of  the  mooring  line. 

Radii  of  the  t u:  ed  spar  buov. 

Spat  1 a 1 vector  in  the  x-y  (horizontal) 
describing  propagation  of  the 
surface  wave. 

i css- s-  ct  i ona  1 area  of  the  surface  buoy. 

C'j  o.  --sect  lonul  areas  of  the  tu'ed  spar 
buoy . 

spatial  coordinate  of  the  mooring  line. 

The  tension  vector. 

.Effective  tension  variable  along  the 
mooring  line. 

Mcner  t of  the  tension  force  f. 

Tension  force  vector  on  a mooring  line  note 
Time  variable. 

Fluid  velocity  vector. 

Relative  fluid  velocity  vector. 

Normal  and  tangential  components  of  UR. 

Ve 1 oc i t y vector. 

Unit  vector  along  the  moorinq  line. 

Array  of  tianslation  and  rotation  rates. 
Velocity  of  water  particle  w in  frame  B. 

Tangential  relative  fluid  velocity  at 
bottom  of  the  tuned  spar. 

Components  of  surface  current. 

Tangential  relative  fluid  velocity  at  step 
of  the  tuned  spar,  also  static  velocity 
of  the  drogued  drifting  buoy  system. 
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1.0  INTRODUCTION 

Oceanographic  measurements  from  moored  buoy 
systems  are  contaminated  by  mooring  motion.  Numerous 
articles  have  been  published  in  the  last  three  years 
pointing  out  quantitatively  the  errors  introduced  by  the 
surface  wave  field.  Large  vertical  excursions  have  been 
experienced  at  great  depths  by  instruments  located  on 
the  synthetic  rope  portions  of  surface-buoyed  mooring 
lines  (WUNSCH  and  DAKLEN,  1974) . Observations  during  the 
POLYMODE  experiment  showed  subsurface  mooring  lines  to 
have  undergone  hundreds  of  meters  of  vertical  excursions. 
For  any  operational  moored  instrument  system  we  need  two 
models  whose  accuracy  can  be  established  with  known 
confidence:  a mooring  system  model  from  which  the  motion 

environment  of  instruments  can  be  predicted,  and  a model 
of  the  instrument  motion  response  characteristics  from 
which  the  measurement  error  can  be  estimated.  Possession 

of  these  models  is  required  for  more  effective  mooring 

\ 

systems  design  and  for  optimum  interpretation  of 

i 

oceanographic  measurements  obtained  from  moored  systems. 
The  first  of  these  two  models  is  the  main  subject  of  I 
this  report. 

An  alternative  to  the  moored  approach  is  the 
use  of  surface-trackable  drogued  drifting  buoys.  Such  a 


buoy  system  employs  a high  drag  device  (or  drogue)  at 
some  depth,  tethered  to  a trackable  buoy  at  the  surface. 

The  major  impediments  in  the  widespread  use  of  this  approach 
have  been  inadequate  component  and  system  design.  Both 
of  these  impediments  can  be  removed  by  the  development  of 
dynamic  modeling  of  drogued  buoys.  A part  of  this  report 
deals  with  the  dynamic  modeling  of  drogued  buoys. 

Dynamic  models  of  drogued  buoys  do  not  exist 
except  in  primitive  form,  while  such  models  of  moored  buoys 
and  free-drifting  buoys  are  further  along.  This  difference 
stems  from  the  fact  that  drogued  buoy  systems  are  dynamically 
very  complex,  and  they  have  only  recently  been  considered 
essential  to  major  programs.  Even  though  many  dynamic 
mathematical  models  for  moored  buoys  and  free-drifting 
buoys  are  available,  not  much  has  been  done  toward  the 
evaluation  of  these  models  using  full-scale  ocean  test 
data.  It  is  the  purpose  of  this  report  to  present  a 
general,  computationally  efficient  approach  for  analyzing 
moored  buoy  systems,  free  drifting  buoys,  and  drogued 
buoy  systems  and  to  simulate  this  analysis  on  the  computer 
so  that  these  models  can  then  be  readily  evaluated  with 
full  scale  ocean  test  data. 

This  report  constitutes  parts  3 and  4 of  a 
four  part  report.  Parts  1 and  2 were  published  in 
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volume  1 and  dealt  with  the  "Three  Dimensional  Static 
Analy. is  and  Design  of  Single  Point  Taut  and  Slack  Moored 
Buoy  Systems."  (CHHABRA,  197  3)  . 

1. 1 Background 

We  at  Charles  Stark  Draper  Laboratory  Inc.  have 
evaluated  one  of  the  mathematical  models  presented  in 
this  report.  The  mathematical  model  of  a subsurface 
mooring  system  (Section  2.2.1  and  simulation  5.2  in  this 
report)  has  been  shown  previously  to  predict  the  mooring 
motion  forced  by  ocean  currents  of  periods  greater  than 
15  min. (CHHABRA,  DAHLEN  and  FROIDEVAUX,  1974).  The  model 
was  evaluated  in  a full  scale  ocean  test  that  provided 
experimental  data  on  mooring  response  and  ocean  current 
forces.  The  ocean  test  was  conducted  jointly  with  Woods 
Hole  Oceanographic  Institution  on  R.  V.  Chain  cruise 
107.  We  obtained  a record  of  the  motion  of  an  acoustic 
transceiver  near  the  top  of  the  subsurface  mooring  line 
at  the  500-m  depth.  The  transceiver  sent  out  a sound 
pulse  every  minute  and  recorded  the  four  return  times 
of  replies  from  four  near-bottom  acoustic  transponders 
at  5460m.  By  comparison  with  this  and  other  data  from 
precision  pressure  recorders,  tensiometers,  and 
inclinometers,  the  mooring  model  had  been  found  to 
predict  well  the  observed  motions.  The  r.m.s.  difference 


between  the  experimental  and  predicted  trajectory  of  the 
acousuic  transceiver  was  11.8m,  about  10*  of  the  mean 
excursion.  This  mathematical  model  was  evaluated  a 
second  time  by  data  from  the  central  mooring  (Mooring 
No.  1,  Station  431)  of  the  Mid-Ocean  Dynamics  Experiment 
(MODE).  In  that  study  (CHHABRA,  1976),  the  current  record 
from  the  topmost  vector-averaging  current  meter  on  the 
central  mooring  of  the  MODE  experiment  was  corrected  for  the 
effects  of  mooring  notion,  and  [ower  spectra  of  the  uncorrected 
and  corrected  signals  were  compared.  The  correction  was 
•:  1 cm  s'1,  for  that  mooring  line  (vertical  excursion 
< 12  m) . Creep  in  the  synthetic  portion  of  the  mooring 
line  was  also  identified. 

In  addition,  we  are  planning  to  evaluate  a 
second  mathematical  model  presented  in  this  report.  The 
mathematical  model  is  of  a tuned  spar  buoy  (35  ft.  long, 

1 ft.  dia.  at  W.L.)  tethered  to  a subsurface  mooring  line 
by  a stiff  buoyant  line.  An  instrument  line  is  hanging 
from  the  base  of  the  spar  (Sections  2.1.1,  2.2.2,  2.3 
and  simulation  5.3).  Thin,  configuration  was  recently 
tested  in  the  ocean  during  the  October  1976  ONR/NDBO 
Mooring  Dynamics  Experiment.  A total  of  fifteen  motion 
sensing  instruments  (4  Force  Vector  Recorders,  6 
Temperature/Pressure  Recorders,  1 POPMIP,  and  4 Acoustic 


Beacons)  were  attached  along  this  mooring  system.  The 
evaluation  and  improvement  of  this  mathematical  model 
would  be  done  by  comparing  measured  responses  with 
those  computed  by  the  computer  simulation  for  the 
measured/observed  environment. 

The  above  mentioned  evaluation  and  improvement 
is  planned  to  be  completed  in  CY77.  If  such  an 
evaluation  and  improvement  is  done,  the  only  mathematical 
model  to  be  evaluated  and  improved  from  the  analysis 
presented  in  this  report  would  be  the  drifting  drogued 
buoys.  We  plan  to  do  that  task  in  the  near  future. 

1.2  Assumptions,  Capabilities,  and  Limitations 

Surface  gravity  waves  treated  in  this  report 
have  a single  frequency  and  amplitude,  propagating  in  a 
single  direction.  It  is  assumed  that  the  surface  buoy 
is  oscillating  in  the  path  of  small  (i.e.  amplitude 
of  wave  train  much  less  than  its  wavelength)  incident 
surface  waves  which  are  long  in  relation  to  the  body 
dimension  in  the  direction  of  wave  propagation.  Wave 
direction  has  been  generalized  to  include  any  arbitrary 
direction  in  the  horizontal  plane.  For  the  three- 
dimensional  analysis  of  surface  buoys,  we  also  _ssume  small 
displacements  of  buoy  when  compared  to  vertical 
dimension  of  the  buoy.  Hence  only  the  first  powers  of 


small  quantities  are  retained  for  the  three-dimersionjtl 
surface  buoy  analysis.  Moments  of  inertia  about  all 
axes  in  a horizontal  plane  are  taken  equal.  Even  though 
moment  of  inertia  about  the  longitudinal  axis  of  the 
spar  buoy  is  negligible,  it  is  included  for  numerical 
computational  purposes. 

In  all  mathematical  models  an  earth  fixed 
frame  is  assumed  to  be  the  valid  inertial  frame.  Viscous 
drag  forces  are  computed  based  on  the  square  drag  law. 

For  cylindrical  shapes  these  forces  are  assumed  to  act 
in  directions  normal  and  tangential  to  the  longitudinal 
axis;  and  for  any  other  shape  they  act  in  the  direction 
of  the  relative  flow.  For  the  case  of  a tuned  spar 
viscous  drag  forces  and  'added  mass'  due  to  the  bottom 
base  and  the  step  are  also  included.  The  ideal  fluid 
damping  (wave-damping)  as  derived  by  Newman  (1963)  for  a 
spar  buoy  was  found  to  be  negligible  and  hence  is  neglected. 

All  mooring  lines  are  considered  elastic. 

In  the  continuous  line  formulation  of  mooring  lines, 

i 

nonlinear  elasticity  dependent  on  the  prior  loading  history 
as  explained  in  volume  1 (CHHABRA,  1973)  is  considered. 

For  the  lumped  parameter  formulation;  a linear  stress-strain 
curve  is  derived  from  the  nonlirjear  curves  for  the  range 
of  stresses  in  study.  Dynamic  effects  on  these  stress- 
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strain  curves  are  neglected.  Internal  damping  forces  in 
the  mooring  lines  are  assumed  negligible  when  compared 
with  viscous  drag  and  stiffness  forces.  The  hydrostatic 
pressure  forces  are  treated  more  rigorously  than  they 
are  in  the  traditional  method  employed  in  most  previous 
studies.  As  it  turned  out  this  new  methou  changed  only 
the  stress  distribution  along  the  mooring  system,  while 
the  configuration  of  the  mooring  system  remained 
essentially  unchanged.  Instruments  or  buoyancy  packages 
attached  along  the  mooring  lines  are  treated  as 
concentrated  forces  which  have  length  and  give  rise  to 
forces  as  explained  in  Section  2.2.  A time  varying 
current  profile  of  any  nature  and  shape  can  be  inputted 
to  the  computer  programs.  A list  of  current  profiles 
used  in  various  simulations  is  given  in  volume  1 
(CHHABRA,  1973)  . In  its  present  form  no  allowance  for 
shrinkage  or  creep  of  the  synthetic  ropes  is  taken  into 
account,  but  creep  was  identified  in  CHHABRA  (1976). 
Elongations  due  to  rotation  of  non-torque  balanced  cables 
is  also  not  considered.  In  the  mooring  line  models 
viscous  drag  forces  due  to  the  velocity  field  generated 
by  the  surface  waves  are  included.  Wave-damping  forces 
are  neglected.  In  the  continuous  line  formulation  of  the 
mooring  line,  exciting  forces  exerted  on  the  line  by 
the  wave  system  are  neglected;  whereas  i.i  the  lumped 


parameter  formulation,  exciting  forces  are  taken  equal 
to  the  *Froude-Krylov'  forces  as  it  is  assumed  that  the 
presence  of  the  mooring  line  does  not  disturb  the  wave 
particle  motion.  For  added  mass  purposes,  it  is  assumed 
for  the  continuous  line  formulation  that  the  body 
(continuous  cylindrical  line)  notion  accelerates  fluid 
only  in  the  direction  normal  to  its  longitudinal  axis. 
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2.0  THEORETICAL  ANALYSIS  - MATHEMATICAL  MODELS 

This  section  derives  all  the  equations  of  motion 
pertinent  to  the  analysis  presented  in  this  report. 

Surface  floats  are  analyzed  in  Subsection  2.1,  mooring  line 
in  2.2,  their  attachment  in  Section  2.3,  and  a window 
shade  drogue  in  Section  2.4.  The  mooring  line  analysis 
includes  all  subsurface  floats,  instruments  etc.  attached 
to  the  mooring  line. 

2 . 1 Surface  Buoys 

As  is  well  known,  the  analysis  of  the  wave 
induced  response  of  floating  Lotties  is,  in  general,  a 
most  formidable  task.  Initially  the  action  of.  the  fluid 
must  be  decomposed  into  real  (viscous)  and  ideal  (inviscid) 
effects.  Each  effect  then  must  be  modeled  as  to  its 
interaction  with  the  floating  body.  Further,  it  is 
convenient  to  assume  that  the  ocean  waves  are  "gentle" 
enough  so  as  to  permit  first  order  linear  surface  wave 
theory  to  seive  as  a foundation  for  the  calculation  of  the 
wave  exciting  forces. 

The  ideal  fluid  problem  is  still  so  difficult 
in  general  that  only  the  simplest  body  shapes  are 
amenable  to  rigorous  analytical  solution  (potential 
flow  theory) . This  solution  involves  the  determination 
of  the  (velocity)  potential  function  4>(S,t)  which  rust 
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satisfy:  1)  Laplace's  equation,  2)  the  kinematic  boundary 

condition  on  the  moving  body  surface,  3)  the  free 
surface  (and  bottom)  boundary  condi tion(s)  and  4)  the 
radiation  condition  at  great  distance  from  the  body. 

Assuming  this  potential,  -1  (53,  t) , has  been  determined 
by  some  means,  the  next  step  is  to  substitute  it  into  the 
Bernoulli's  expression  for  fluid  pressure  which  in  turn 
is  integrated  over  the  immersed  portion  of  the  body  to 
yield  the  instantaneous  force  and  moment  vectors.  In 
principle  this  is  a straight-'orward  procedure.  In  practice 
it  is  hardly  ever  possible  to  find  a tractable  potential 
function  which  satisfies  the  above  four  conditions.  In  a 
few  notable  cases,  however,  for  simple  geometries  and 
small  body  motions,  rigorous  solutions  have  been  worked 
out.  In  particular  Newman  (1963)  has  derived  the  linearized 
equations  of  motion  for  a vertical,  cylindrical  * ar  buoy 
responding  to  the  influence  of  a unidirectional  wave 
train.  The  ax i symmetry  of  the  buoy  as  well  as  its 
postulated  slenderness  was  greatly  exploited  in  the  work 
to  yield  managablo  results.  In  Newman's  work  as  well  as 
in  others  of  comparable  rigor,  the  computed  potential 
is  the  result  of  an  intricate  distribution  of  singularities 
(sources,  sinks,  dipoles,  etc.)  within  or  over  the  wetted 
body  surface. 

From  these  studies,  it  turns  out  that  the  linearized 
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hydrodynamic  forces  and  moments  in  general  body  geometries 
may  be  decomposed  into  constituents  proportional  to 
body  acceleration,  body  velocity,  displaced  fluid 
acceleration,  and  the  displaced  fluid  velocity.  In 
general,  rigorous  integration  of  the  fluid  pressure,  both 
hydrostatic  and  wave,  over  the  body  surface  yields: 

1)  the  'Froude-Krylov"  force  from  the  pressure  distribution 
due  to  the  undisturbed  wave  system;  2)  the  "diffraction" 
force  from  disturbance  of  the  waves  by  the  presence  of 
the  body;  3)  the  force  due  to  the  motion  of  the  body; 
and  4)  the  hydrostatic  pressure  force.  The  Froude-Krylov 
force  and  the  diffraction  force  in  combination  are  also 
known  as  the  exciting  forces  exerted  on  the  body  by  the 
wave  system.  The  Froude-Krylov  force  equals  in  value 
to  the  product  of  the  mass  of  the  displaced  fluid  times 
the  acceleration  cf  the  local  undisturbed  fluid  particles. 
The  diffraction  force  and  the  force  due  to  the  motion 
of  the  body  both  yield  "added  mass"  coefficients 
proportional  to  accelerations  and  wave-damping 
coefficients  proportional  to  velocities  of  the  fluid 
and  the  body  respectively.  As  shown  in  Chung  (1976) , 
the  added  mass  coefficients  are  same  (opposite  signs)  for 
both  the  fluid  and  body  accelerations.  Also  the  two  wave- 
damping coefficients  have  the  same  (opposite  signs)  value. 
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Hence,  these  two  forces  may  be  combined  to  represent 
an  "added  mass"  force  proportional  to  the  relative 
acceleration  (fluid  acceleration  minus  the  body  acceleration), 
and  a wave-damping  force  proportional  to  the  relative 
velocity.  Other  forces  (apart  from,  the  fluid  pressure 
forces)  acting  on  the  body  are:  1)  the  weight  of  the  body; 

and  2)  the  actual  unbalanced  force  which  accelerates  the 
body.  In  the  real  (viscous)  fluid  another  force,  the 
viscous  drag  force,  acts  on  the  submerged  portion  of  the 
body . 

Analysis  in  this  section  includes  all  the  above 
mentioned  forces,  as  derived  from  the  potential  flow  theory, 
except  the  following  deviations:  1)  the  ideal  fluid 

damping  (wave-damping)  is  omitted;  2)  the  added  mass 
force  is  treated  slightly  differently,  and  3)  the  viscous 
drag  forces  are  added  on  to  the  general  equations  of  motion. 

The  ideal  fluid  damping  is  omitted  here  for 
expediency.  Later  analysis  of  specific  problem  addressed 
in  this  report  showed  the  wave  damping  force,  as  derived 
by  Newman  (1963)  for  a spar  buoy,  to  be  negligible 
compared  to  the  viscous  drag  force  on  the  submerged  body. 

The  added  mass  force,  which  is  proportional  to  the 
relative  acceleration  is  separated  into  two  terms;  an 
added  mass  force  proportional  to  the  longitudinal  component 
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of  the  lelat-.ive  acceleration,  and  an  added  mass  force 
proportional  to  the  transverse  component  of  the  relative 
acceleration.  Each  of  these  two  terms  is  multiplied  by 
a different  constant  coefficient,  depending  on  the  shape 
of  the  floating  body.  These  coefficients  and  the 
co'-'ff icients  used  to  derive  the  viscous  drag  force  are 
experimentally  determined  hydrodynamic  coefficients. 

As  will  be  shown  below,  our  general  equations  of 
motion  reduce  to  the  ones  given  in  Rudnick  (1967)  for  a 
particular  value  of  these  coefficients.  Rudnick  in  his 
analysis  used  the  reasoning  of  Lamb  (1945)  who  treats 
a uniform  two-dimensional  flow  across  a long  circular 
cylinder. 


Formulation  of  the  General  Equations  of  Motion: 
Derivation  of  the  equations  of  motion  for  surface  piercing 
buoys  in  a train  of  regular  harmonic  ocean  waves  is 
given.  In  addition  a surface  current  (not  due  to  surface 
waves)  is  present.  Surface  buoys  are  considered  as  rigid 
bodies  with  six  degrees  of  freedom.  The  problem  under 
consideration  is  represented  schematically  in  Figure  2.1. 
An  earth  fixed  car tesian coordinate  system  (z  positive 
upwards)  is  situated  at  the  undisturbed  level  of  the  free 
surface.  Call  it  the  frame  E with  origin  at  o; 
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Figure  2.1  Coordinate  System  for  Surface 
Piercing  Buoys 


A /\  A 

and  and  iz  the  unit  vectors  along  x,y,  and  z 

respectively.  Frame  B is  a body  (surface  buoy)  fixed 
cartesian  coordinate  system  with  its  origin  at  the  centroid 
(c)  of  the  surface  buoy,  and  x',y',z'  being  parallel  to 
x,y,z  respectively  when  the  ouoy  is  in  the  upright  and 
non-rotating  position.  In  this  analysis  R represents 
a displacement  vector;  v(  a velocity  vector,  and  A,  an 
acceleration  vector.  It  is  assumed  that  the  surface  buoy 
is  oscillating  in  the  path  of  small  (i.e.,  amplitude 
of  wave  train  much  less  than  its  wavelength)  incident 
surface  waves  which  are  long  in  relation  to  the  body 
dimensions  in  the  direction  of  wave  propagation.  Point 
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p is  the  location  of  the  centroid  of  a differential 
disk  (height  = dz')  of  the  surface  buoy  (Figure  2.1). 

First,  the  wave  direction  will  be  generalized 
to  include  any  arbitrary  direction  in  the  x-y  plane.  To 
this  end  let  the  wavenumber  vector  K specify  the  direction 
of  propagation  (Figure  2.2)  where: 

_ U>‘  A A 

K = — (cosp  iv  + sinS  i„) 

g y 

Z 


(wave  direction) 


Figure  2.2  - Generalized  Wave  Direction 

The  generalized  velocity  potential  field  of  the 
incident  wave  system  may  now  be  written: 

I v | 7 

■5  = — - e 1 * 1 z cos  (K.S  - cot) 

CO 

A A 

where;  s’  = i x + i y;  and  the  wave  surface  elevation  £ is 
x y 
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given  by: 


C = - sin  (K.f  - uit) 


Here,  g is  the  acceleration  due  to  gravity,  and 
is  the  amplitude  of  the  incident  wave  of  frequency  u. 
In  the  earth  fixed  frame  (frame  E) , assumed  to  be  a valid 
inertial  frame  for  this  problem,  we  define  the  fluid 
particle  acceleration  vector: 


A =i  <5  + i <5  + i ‘j> 

Ew  x vx  y vy  * : 


where  w is  the  particle  of  water  next  to  the  differential 
disk  of  centroid  p.  It  is  assmed  that  A^  is  constant 
over  the  entire  differential  disk.  Components  of  A£w  at 
point  p are  given  by: 


d 3$ 

dt  3x  ~ g,'x^o 


cos (K. C 


<i>t) 


_d_  _3$ 
dt  3y 


K|zw 


y o 


cos  (IC.T 


cot ) 


and  *z  = H = g!Kl  ’O  e^K^Zw  sin  (K.E-  - ut) 

where;  Kx  = |k|  cos3;  Ky  = | K j sinB,  and  the  higher  order 
terms  in  $2  have  been  neglected. 
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IIow;  R = R + R 

op  oc  cp 


% ■ tRopJE  - t«oclE  * WEB  * Rcp 


and 


AEp  = ^Roc^E  + ^-^EB^E  x Rcp  + WEB  X (WEB  X Rcp* 


where;  WEB  is  the  angular  velocity  of  the  buoy. 

Let  us  also  define  the  relative  acceleration  vector 
of  the  fluid  particles  with  respect  to  the  buoy  (AR) , at 
the  location  of  point  p as: 


AR  ~ AEw  “ AEp 


A A 


C-.A80  let,  Art  = (Ar  . i z i ) i z » 


and 


ARN  Ar  ~ art 


The  vector  force  equations  of  motion  may  now 
be  written  as: 


f"1 

f“2 

dn  AEp  = 
i 

dB  AEw  + CN 

£ 

(U2 

u3 

CT 

dB  Art  - 

dB  g 

£ 

u. 


dB  A^  + 


dM  g 


(2.1) 
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wliere  the  viscous  drag  forces  and  the  tension  forces  (due 
to  the  attached  mooring  lines)  are  left  out  for  later 
introduction. 

Here;  £ = " zc,  U1  = H - ZQ,  u2  = hQ  - 7 Q,  and 
u^  = hi  - 7.Q.  7-c  is  the  distance  between  c and  the  bottom 
of  the  buoy.  H is  the  length  of  the  buoy,  hQ  is  the 
instantaneous  draft  of  the  buoy  measured  up  to  the  mean 
free  surface,  and  h.^  is  hQ  plus  the  wave  elevation  component. 

A 

Distances  7.  , SI,  h„  and  h.  are  measured  along  i dM 
is  the  mass  of  the  differential  disk  (Figure  2.1)  of 
height  dz'  and  dB  is  the  mass  of  water  displaced  by  this 
disk.  CN  AND  CT  are  the  appropriate  hydrodynamic  constants. 

In  equation  (2.1)  the  integrals  in  order  of  their 
appearance  represent;  1)  the  actual  unbalanced  force  which 
accelerates  the  buoy;  2)  the  Froude-Krylov  force; 

3)  the  normal  component  of  the  added  mass  force;  4)  the 
tangential  component  of  the  added  mass  force;  3)  the 
hydrostatic  pressure  force;  and  6)  the  weight  or  the 
buoy.  This  equation  reduces  to  the  equation  in  Rudnick  (1967), 
which  was  derived  for  a spar  buoy,  for  C„  = 1.0,  and 

cT  = 0.0. 

The  vector  moment  equation  can  now  be  written 
similarly  as:  (moments  about  C.G.) 


/ 


/ 


\ 

/ 


/ 
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dM(Rcp  X ACp) 


dn(Rcp  x AEW)  + 


+ cn  JdB(Rcp  x aki)  + 


+ CT  jdB(Rcp  X Art) 


(2.2) 


- jdB(Rcp  x g)  + |c:m(Rcp  x g) 


It  is  also  necessary  to  determine  an  appropriate 
transformation  scheme  to  relate  vectors  in  the  E frame 
to  vectors  in  the  B frame.  A sequence  of  rotations 
(Figure  2. 3)  about  the  body  x',y',z'  axes  (roll  = 
pitch  = n , ar.d  yaw  = y)  respectively  yield  the  following 
transformation . 


— 

— 

A 1 

i , 

X 

cos '3  cosy 

sinF  cos: 

sin?  sin? 

ix 

+ cost1  sinO  sin? 

-cos ? sinO 

COSi $ 

i » 

V 1 

-sinC'  cosO 

cosy  cos? 

cosi3  sino 

a 

iv 

y 

-sinO  sinO  sin? 

+sin?  sinO 

C 050 

y 

v. 

sin^' 

-cos9  sin? 

COS  6 COS? 

iz 

— 

The  inverse  of  this  transformation  matrix  is  given 
by  its  transpose.  The  force  and  moment  equations  (2.1) 
and  (2.2)  may  be  combined  and  written  in  matrix  form  as: 


Step  2 - Pitch  (0) : Rotate  x 


Step  3 - Yaw  ( |>)  : Rotate  X£#y 


,2^  about  to  obtain  X2»y 


about  Z2  to  obtain  x',y'. 


Figure  2.3  Transformation  Matrix  Rotations 
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[Mj  • gj  IV,  - {Fnt  + (F0:  ♦ IFft) 


(2.3) 


Here  the  matrix  [M]  is  composed  of  the  tensors 
of  inertia  and  added  inertia.  The  f V } column  vector  has 
three  translation  rates  and  three  angular  rates.  The  vector 
(FDi  contains  the  wave  exciting  forces  and  moments.  The 
vector  f Fq } comprises  the  hydrostatic  pressure  and  weight 
restoring  forces  and  moments  (i.e.  gravitational  and 
buoyancy  effects);  and  {F^}  is  any  additional  forces 
and  moments  including  viscous  drag  forces  and  moments 
introduced  next  and  tension  forces/moments  from  attached 
mooring  lines  introduced  in  Section  2.3. 


Viscous  drag  forces  and  moments  • A square 
drag  law,  where  viscous  drag  forces  are  proportional  to 
the  square  of  the  relative  fluid  veloicty  is  used.  In 
Figure  2.1 


R — R + R + P 
ow  oc  cp  pw 


VEw  - VEc  * WEB  * Rcp  + lRp„]E 


If  point  p and  w are  overlapping  then; 


^Rpw^E  ~ ^Rpw^B  VBw  VEw  VEc  ~ WEB  X Rcp 


-??  - 


Here;  V 


Ew 


i (V  + i)  + i (V  + $ ) + 
x ox  x y y y 


i (V  + <\  ) and  V_. , is  the  relative  velocity  of  water 

2 O Z Z BW 

seen  by  the  buoy  at  the  local  point  p.  'Hence  tho 
viscous  drag  force  on  the  submerged  portion  of  the  buoy  is 
given  by: 


DF 


- * ^ i: 


dAD  i VBw  1 1 <1\: 


(2.4) 


where,  p is  the  water  density,  CQ  the  drag  coefficient,  and 
dAQ  the  appropriate  differential  area.  The  viscous  drag 
moment  can  similarly  be  written  as: 


DM  = 


I 0 CD 


dVVBw!  <Rcp  * VB„> 


(2.5) 


Limits  of  integration  for  both  (2.4)  and  (2.5) 
are  the  entire  submerged  depth  of  the  buoy. 

In  the  remainder  of  Section  2.1  the  equations 
(2.1),  (2»2)>-  (2.4)  , and  (2.5)  presented  above  will  be 
specialized  to  the  specific  problems  at  hand.  The 
matrices  of  equation  (2.3)  will  then  be  derived  and 
presented  for  these  specific  problems. 


2.1.1  Spar  Buoy 

A tuned  spar  buoy  (Figure  2.4)  will  be  analyzed 
in  this  section  to  obtain  general  spar  buoy  equations. 


■ I 

I 


/ 


/ 


! 


Cylindrical  spar  buoy  equations  car  then  be  obtained  as  a 
special  case  of  the  tuned  spar  buny . Doth  two-dimensional 
and:  three-din*  ns lona  1 analysis  will  bo  presented.  *We 
will  find  the  equations  irnorion  of  its  response  for  a 
tram  cf  surface  gravity  waves  having  a single  frequency 
arc  amplitude,  propagating  m any  s i g le  direction. 
Variation  in  wave  acceleration  ar.d  velocity  over  the 
horizontal  buoy  dimensions  are  neglected.  In  addition  in 
the  three-dimensional  analysis,  all  rotational  angles 
(roll,  pitch,  and  yaw)  arc  restricted  to  be  small. 

Two-Dimensional  Analysis:  In  this  analysis  K 

is  a scaler  with  :■  = 0;  and  fl  is  replaced  by  the  scaler 

x which  is  the  direction  of  propagation.  A so,  X^w  = 

. • 
l ; + i ...  . 

XX  z z 


Vi  ■ IAr 


V 1 V 


and  the  appropriate  transformation  with  i - 0,  and  < = 0 

becomes 


I ix' 


. > - 


' r 

; cos  9 -sinO 

l 

Ax 

- 1 

;sinrt  cosh* 

i | 

17.  | 

\ i 

-2L- 


Now  let. 


and 


Also, 


Define; 


and 


Integrals 


1. 


A • « 


'■ Soc]E  - xc  + ^ "c 


wcb  - y 


Cp  ~ p Z* 


__  ■ ••  • 2 

A£p  = ix  [:<c  -»  Zp  (0  cost' -6  sin1)] 

~ ..  • 2 

+ i_  [y.  + z'  (-0  sin0-9  cosO)] 

Z C p 


A£w  = 2fo  e*Zv  ^ ix  cos(Kxc-  + i2  sir.{Kxc  - t)  ] 


I U2 


Pi  " 


Qi  = 


ru2 


’ * 0,1,2 


dB(2pX 

Kzw  i 

dBc  (z^r  i = 0,1 


of  equation  (2.1)  can  now  be  written,  as; 


fdM  A.,  - M[i  x + i z ] 

J tp  x c z c 


Where  M is  the  total  mass  of  the  buoy. 


dB  Anw 


= ' '..QoUj  coslK*c  - ■- 1)  . i 2 oin(Kxc  - .t)l 
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3«  CN  \ dB  ARN  = CN  ix^’^0Qo[cOS(Kxc  “ ')t:)  COs20  ’ 

- sin(Kxc  - ->t ) sinO  cosO]  + 

••  •*  2 
+ p (z  sinO  cosO  - >:  cos  0)  - Pi  cosO 9}  + 

o c c J- 

2 2 

+ C„  i {(jj  C„Q„[sin  (Kx„  - <.,t)sin  0 - cos  (Kx  - t) 

N 2 O O C C* 

2 

sinO  cosQ]  + (::  sinO  cosO  - z sin  0 ) + 

o c ^ 

+ P^sinOO} 

4.  CT  dB  Art  = CT  ix{u)2rJOQ0Ccos  (Kxc  - it)  sin29  + 

+ sin(Kx  - <’t)sin6  cosO]  + 

v 

••  2 •• 

+ PQ(-xc  sin  G - 7-c  sinO  cosft)  + 

^ A « 

+ P^sinOGi  } + CT  iz(w  C0Q0[cos(Kxc  - t)sin0cosr  + 

+ sin(Kxc  - u)t)cos20]  + PQ  {-xcsin6cos0  - 

**  2 * 2 

- B cos^G)  + P^  cosOO  } 


5.  - 

[u3 

dB  g = g 

rU2  . 

dB  i2  + g 

->;Qsin  (Kxc~..>t) /cosO 
dB  iz 

. 

l 

l 

o 

Also;  dB  = pSQdz',  where  SQ  is  the  cross-sectional 
area  of  the  disk  of  height  dz'.  Then; 


i 
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f u- 


dB  g = iz[Pog  ~ cS02gr-0  sin(Kxc  ~ «*’t)/cos8] 


where;  So2  - v ^ 


6. 


dMg  = - Mg  i . 


Integrals  of  equation  (2.2)  can  also  be  written 


as: 


dM(Rcp  x V - [M  2'2  e iy  - V 6 iy 


— — * 2 

2.  dB(R_.p  x aew>  = iy^  r'0Q].fcos  (Kxc-u)t)cos0  - 


- sin(Kx  -.,t)sin9]} 

V 


3.  C 


N 


dB(Rcp  x ARN5  = CN  iy{<i)  ^oQi[cos(Kxc-;Jt)cos0 
- sin (Kxc- . t) sinP ] - P^cosO  xc  + 


+ P,sinO  7.  - P,0) 

1 C 2 


4.  CT  dB(R  x Art)  = 0 


fU3 


5.  - 


ru- 


dB(Rcp  x g)  = - 


dB(RCp  x g) 


= iyt-Piy  sine) 
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6. 


dM(Rcp  x g)  - 


0 


Viscous  Drag  Forces  and  Moments:  /or  a tuned 

opar  buoy  the  viscous  drag  forces  are  assumed  to  be 
acting  in  normal  and  tangential  directions  proportional  to 
their  respective  drag  coefficients  and  areas. 

V ' Mvox  + ■ *c  - S6co,9]  + 

A 0 , 

+ iZ^OZ  + $Z  “ ZC  + Z^,0Sin0] 

= *x'C<Vox  + *x  - *C)COS0 

- (v02  + *z  - zc)sin0  - Z^e] 

+ V^ox  + “ *c)sin0 

‘ + {Voz  + *z  “ zc)cos0] 

“ ix'  VBwx ' + iz'  VBwz* 

Drag  force  components  can  now  be  written: 


P^"DN j r I VBwx ' I VBwx‘ 


pCDT71 


rl VBwz' lVBwz' 


dz ' 
dz ' 


* 7 ‘1V{rl|vb|vb  + 'ri  - 4>lvs|vs) 


Here;  CDN  and  CDT  are  normal  and  tangential  drag 
coefficients.  CDP  is  the  pressure  drag  coefficient  due 
to  the  bottom,  and  the  step,  of  the  tuned  spar  buoy. 


vs  - Vz'  (z'=  hl  - =c> 


and 


v = V . ( * ' = -1  ) 
vb  vBwz  v c‘ 


Also;  DM  = pC 


DN 


rlVBwx'lVBwx’  Zp  dz’ 


Using  the  transformation  matrix: 


DF. 


DF, 


cos9  sinO 


-sinQ  cos0 


DF  *• 


DF  z, 


Added  Mass  Force  Due  to  the  Step  and  Bottom  Base 

of  the  Tuned  Spar:  In  addition  to  the  forces  mentioned 

above,  an  added  mass  force  proportional  to  the  tangential 

body  acceleration  (AEt)  and  acting  perpendicular  to  the 

step  and  the  bottom  end  of  the  tuned  spar  is  considered. 

Lamb  (1945)  in  his  analysis  of  a cylindrical  spar  moving 

in  still  water,  gives  an  added  mass  coefficient  equal 
3 

to  4/3pr  . Following  representation  is  used  in  this 
analysis  for  this  force  (AM) . 


-3C  - 


AM  - j ap 


rl  AEtl  2 ’=-7, 


(rl-r2>3  AEtlz’=h,-r 

1 c 


(2.6) 


where;  AEt  = (AEp* ig , ) iz , , and  a is  a constant. 

AM  is  added  to  the  left  hand  side  of  equation  (2.1) 

R „ x AM  would  be  equal  to  zero, 
cp 

Combining  all  these  forces  and  moments;  the  matrix 
coefficients  of  equation  (2.3)  can  be  written  as: 


2 2 

= M + Pq(Cncos  6 + CTsin  0)  + 

4 «pCr13+  (ri_r 2)  3]sin2e 

Ml2  = PQsin0cos0 (CT  - CN)  + 

4 3 3 

apLr^t  (r1-r2)  ]sin6cosfl 


M13  " CNPlCos6 


M21  = Posin0cos9  (CT  - CfJ)  + 


4 3 3 

uf.[r^  + (r,-r2)  ]sin9cos9 


M 


22 


= M + PQ(CNsin29  + CTcos29) 

3 ap[r13  + (r  -r  ) ^cos*^ 

1 2 


+ 


- / 


-31- 


M 


23 


- C„P,sin0 
N 1 


M 


31 


cnp ^cose 


M 


32 


- C,,P.sin0 
N 1 


H33  " Jyy  + CNP2 


{V}  = Lx_,  7.  , 6];  which  is  a column  vector. 

v C 


2 2 2 
FD1  = W C0Q0tcos  (Kxc-uit)  {l+CNcos  0 + CTsin  0)  + 


sin(Kxc~ut)  sin0cos0 (CT-CN) ] 


2 2 2 
Fd2  = W C0QQ[sin  (Kxc-(.)t)  ( l+CNsin  0+CTcos  6)  + 


cos  (Kx  -ot)  sin0cos0  (CT~CN)  ] 


/ 


Fd3  = u CoQ1[cos(Kxc-wt)cos0-sin(Kxc-',t)sin0]  (1+CN) 


1 Gl  = °* 


\ 

V 


Fg2  = Pog  - P so2  gCosin(Kxc-i..t}/cos0-Mg 


Fg3  = - P^g  sin0 
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O <»  a ^ o 

FA1  = CTP^sinO 0 + DFX  - -j  apt^  +(r1-r2)  ]7.csin9e 

FA2  = C^jCOseo  + DFz  + ^ ap[rx  +(r1-r2)J] 
(h1-2c)cos062 

and  Fa2  = DM 


Three-Dimensional  Analysis;  For  small  angles  the 
appropriate  transformation  matrix  reduces  to: 


A 

V 

— 

1 v -e 

’ A 

*x 

A 

i 

J 

sr 
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A 

1y 

A 

e -<t>  l 

A 

define: 


In  this  analysis  we  let  CN  = 1.0,  CT  = 0.0,  and 


^Roc^E  ^xf'c  + *y  yc  + is'  c 


WEB  * V * V * lz* 


and 


Rcp  zp  ■*’z' 


. . At.„  = 1 


l..O-  * * iy'-'c  - * Vc 


Ep  X c p 


and; 


AEw  oe 


i r.  i .. 

2.  'K-w 


[ i cosr cos (K • f-  t)  + 


as : 


i sin. cos  (K • f-  t)  + i sin(K-r-.t)] 

y z 


Integrals  of  equation  (2.1)  can  now  be  written 


1. 


(iM  A = M ( l + i V + i„z_ 
JTp  x c yJc  z c 


2. 


dB  A, 


Ew 


f Q [i  cosrcos  (K*f>  t)  + 
O O X 


i sin  cos  (K • B-  t)  + i,  sin(K‘f-  t)] 

y z 


also;  Arn  ^x^'x  ” "c  ’ ! ' 


+ + izl"iic  " ;Vc) 


y y c 


3.  jdB  Arn  = ixl  *i2-0cos:-:Q0cos(K«r-  t)  - PD*C  " P]/ 


ro"Zc  ' + £y  a * 0S 1 n ' 0oc 03  ( K • S-  . t ) 


4. 


P y +■  ;» 
o c 1 


_ 3 V o j 

- I dB  g = j dB  cj  + JdB  g i 


P •:  z ) + i P ( <x 

O C 7 O'  C 


i v i 
* c 


sin  ( K • S - -t) 
o 


= iz|Poq  " (So2qfosin  <K*S  - -t)J 


-V- 


5. 


dM  g = - Mg  i. 


Similarly  integrals  of  equation  (2.2)  can  be 
written  as: 


1. 


x V = ix  ^x*  + Ay  V 


where;  1^  » Iyy 


dM  s?'* 


2.  jdB(Rcp  :c  A^)  =-  ix{-u)2f,0Q1sinpcos(K*C-c  t)  } 

+ iy  (w2CcQ1cosPcos  (K*S-  't)> 


dB(Rcp  x = ix { -u)  f^Q^inPcos  (K*  "-0t) +P17c-P24> 


+ P]_#zc}  + iy{u  CqQjCOsBcos  (K*S-.^*:) 


" pixc  - P2°  + Pl0zc}  ' iZPl{0i'c  + *xc} 


u- 


dB(RCp  x g)  = - 


u. 


dB(Rcp  x 9> 


ix(-Pl<>g)  + i (-P^6g) 


JdM(Rcp  x g)  = 


5 


0 


Viscous  drag  Forces  and  Moments : 


Hence; 


and. 


V = i I V + $ - - c ' 0 ] + 

Bw  xL  ox  x c p 

A • • > 

i [y  + <J>  - y + z • <+)  J + 

yL  °y  y c p 


i [v  + ♦ - z ] 

Z L 02  Z CJ 


- V tVox+^x”*c”Zp9+Voy^~^c^~Voz9+zc9^ 

+ iy  • t“vox*+,;,iC+voy+*y’Vzp*+voz*-*C^3 

+ iz  • ^ox^'V-V^^c^V^-y 

/\  A A 

“ Ax'  VBwx*  + iy'  VBwv ' + iz’  VBwz' 


DFx’  = pCDM | - J'Bwx 


r ( VBwx'  + VBwy'  ^VBwx'^dZ 


DFy,  " °CDN  VBWX ' 2 + VBwv'2  (VBwv,)dz' 


’ Bwy 1 Bwy ' 


t>P2.  - f-cDT. 


rlVBW2’lVBwz'  dj>  * 

I eCDp”lrl2lvblvb  * <rl2  - r22'|Vs|Vs) 


Using  the  transformation  matrix: 


26- 


Also 


Therefore: 


[DFx 

1 -v 

e 

DFx' 

-.A-.,,.,- 

O 

*< 

• = 

1 

-<■ 

DFy'  ■ 

iDfy 

-e  <p 

1 

DF  ' 
2 

DM  = 

/«% 

xdDF) 

V 

(DMx.)  + 

V 

(DMyt) 

DMx 

1 - t 

9 

DM  , 
x' 

DMy 

a 

4>  1 

-<P 

DM 

y 

DMZ 

-9  4> 

1 

0 

Combining  all  these  forces  and  equation  (2.6)  of  additional 
added  mass,  we  can  write  the  non-zero  elements  of  matrices 
in  equation  (2.3)  as: 


Mii  = 

M + p0  = M22 

M13  = 

4/3  ap[r13  + (r1-r2)3]9  - 

Pc9 

M15  = 

P1  = M51 

M23  = 

P04>  ~ 4/3  ap [r^3  + (rj-r2) 

i3]^ 

M24  = 

- P1  " M42 

= M 


31 


= M + 4/3  ap[r-j^  + (r^-^)^] 


M43  =~  Pl^ 


M44  = Txx  + P2 


M53  ='  Pl9 


M55  " Iyy  + P2 


M61  = Pl* 


M62  = Pl6 


{V}  = [::  , y , z , 0,  tj>]  is  a column  vector, 

w c c 


FD1  = 2a)  5oQoC0s6cos  (K“S-a)t) 


Fd2  ® 2a)  CoQosin0cos  (K*S-ut) 


FD3  = w S„Q„sin  (K*  R-ait) 


’o“o 


Fd4  «-  2u  ^QQ-LSinBcos  {K«S-a)t} 


FD5  = 2a)  cosBcos  (K*r-wt) 


’owl 


\> 

\ 


FG3  - Po9  - •'  So2gf.0sin(K-C-.„t)  - Hg 

FG4  " - PP9 
FG5  ' " Pl"9 

and  IF.}  = [DF  , DF  ,,  DF  . DM  , DM,.  DM  ] is  a column 

r\  X y Z X y Z 

vector. 

As  can  be  seer  from  these  matrix  elements,  the 
sixth  degree  of  freedom  corresponding  to  ^ does  not  drop 
out;  but  equals  zero.  Mgg  is  introduced  in  the  equations 

for  computational  purposes.  f’*6g  is  given  by: 

M66  = T^rl  hl  + r2 

vrhich  represents  the  moment  of  inertia  of  a tuned  spar 
about  its  longitudinal  axis. 

2.1.2  Other  Shapes 

To  obtain  the  velocity  potential  function  : 
for  any  shapes  other  than  for  simple  geometries  and  small 
body  motions  is  hardly  ever  possible.  lienee  to  obtain 
any  reasonable  solution  to  this  problem  one  has  to  depend 
on  emperical  approaches.  Solution  to  these  problems  has 
thus  been  left  out  of  this  report.  In  order  to  use  the 


i — r 
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analysis  and  computer  programs  of  this  report,  the  readers 
will  have  to  substitute  their  own  elements  for  [m],  (F^), 

{ F^.} , and  {F^}  matrices  of  equation  2.3;  pertinent  to  the 
particular  surface  buoy  in  question. 

2.2  Mooring  Line 

A mooring  line  connects  a surface  or  a subsurface 
buoy  to  the  anchor.  In  this  report,  a line  attached  to  a 
buoy  but  not  to  an  anchor  is  also  considered  a mooring  line. 
Such  a line  could  be  connecting  a drogue  with  the  surface 
buoy  or  be  an  instrument  line  hanging  from  a moored  surface 
buoy.  In  general  a mooring  line  is  made  of  any  type  or 
number  of  materials  (steel,  nylon,  dacron,  etc.)  and  has 
any  type  or  number  of  instruments  (including  subsurface 
floats)  inserted  along  its  length. 

The  mathematical  model  of  mooring  line  dynamics 
will  be  formulated  in  two  different  approaches.  The  "first" 
approach  is  called  the  "continuous  line  formulation".  In 
this  approach  the  mooring  line  differential  equations, 
with  respect  to  the  spatial  coordinate  s,  are  integrated 
incrementally  down  the  mooring  line,  to  obtain  its 
dynamic  equilibrium  at  any  instant  of  time.  Velocities 
and  accelerations  of  the  mooring  line  differential  elements 
and  the  instruments  (including  subsurface  floats)  inserted 
in  the  mooring  line  are  computed  by  differentiating 


\ 


■ 4 . 

i ' 


i\ 

, ■ \ 
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positions found  by  the  dynamic  equilibriums.  This  approach 
was  used  to  model  the  low  frequency  motion  of  a subsurface 
mooring  system,  as  presented  later  in  the  report.  Exciting 
forces  exerted  on  the  mooring  line  by  the  wave  system  are 
neglected  and  so  are  the  wave-damping  forces.  Viscous 
drag  effects  due  to  the  velocity  field  generated  by  the 
surface  waves  is  included.  This  approach  was  found  to  be 
computationally  inefficient,  and  hard  to  solve  numerically 
(due  to  the  differentiation  of  positions  to  find  velocities 
and  accelerations)  for  the  high  frequency  motion  of  a 
surface  moored  system. 

For  the  high  frequency  motion  of  a surface  moored 
system  a "second"  approach  called  the  "lumped  parc^ter 
formulation"  is  presented.  Here  exciting  forces  exerted 
on  the  mooring  line  by  the  wave  system  are  taken  equal  to  the 
"Froude-Krylov"  forces,  as  it  is  assumed  that  the  presence 
of  the  mooring  line  does  not  disturb  the  wave  particle 
motion.  Aqain,  viscous  draa  effects  due  to  the  velocity 
field  generated  by  the  surface  waves  are  included,  and  the 
wave-damping  forces  are  neglected.  1.  both  approaches,  a 
velocity  profile  (could-be  time  varying)  can  be  present 
along  with  the  surface  wave.  Both  formulations  are  presented 
in  three-dimensions  and  can  be  reduced  to  two-dimensions 
when  needed. 
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2.2.1  Continuous  Line  Formulation 

On  a differential  element  of  a continuous  mooring 
line  the  forces  acting  are:  (a)  the  constant  force  due 

to  gravitational  attraction,  (b 1 the  variable  tensile 
forces  transmitted  from  uhe  adjoining  elements,  and 
( c ) the  variable  pressure  (normal)  and  shear  (tangential) 
forces  applied  by  the  fluid.  The  fluid  forces  can  be 
broken  down  into  (1)  the  hydrostatic  pressure  force, 

(2)  the  pressure  force  due  to  the  acceleration  of  the  fluid 
by  the  element  (the  so-called  added  mass  force) , and 

(3)  the  viscous  drag  forces  due  to  fluid  relative  velocity, 
which  have  both  pressure  (normal)  and  shear  (tangential) 
components.  Internal  damping  forces  are  neglected  in 
this  analysis  as  these  are  assumed  small  compared  to 
tensile  and  viscous  drag  forces.  Exciting  and  damping 
forces  due  to  the  wave  system  are  also  neglected.  From 
Newton's  second  law  of  mechanics  the  above  forces  should 
equal  the  mass  of  the  differential  element  multiplied 

by  its  acceleration.  By  computing  these  forces  the 
differential  equations  for  the  mooring  line  dyramic 
equilibrium  are  derived  as  follows: 

The  mooring  line  is  considered  to  be  a cylindrical 
slender  body.  A free  body  diagram  of  a differential 
element  of  length  ds  of  the  mooring  line  is  shown  in 
Figure  2.5. 
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Figure  2.5  Mooring  Line  Differential 
Element 


The  internal  tension  and  the  inclination  of  the 
line  segment  change  to  keep  all  the  above-mentioned  forces 

A 

in  equilibrium.  In  Figure  2.5  v is  a unit  vector  along 
the  mooring  line  given  by: 


/V  A 

= ixcos<|.1  + iyCOS$2  + 


V 
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7 is  a tension  vector,  and  U is  a fluid  velocity  vector 
(obtained  from  the  current  profile  and  the  velocity  field 
generated  by  the  surface  wave) . The  relative  velocity 
vector  UR  is  given  by: 

ur  = u - ? 

c 

AAA 

where;  U = ixUx  + iyUy  + izU2 
.1  * . A . * . 

and,  Y = i x + i,,y_  + i,2_  is  the  velocity  of  point 

c x c y c z c 

c on  the  differential  segment. 

The  viscous  drag  forces  are  computed  according 
to  the  square  drag  law  and  are  assumed  to  act  in  normal 
and  tangential  directions  to  the  cylindrical  line, 
proportional  to  their  respective  drag  coefficients  and 
areas. 

Let:  UR  = UTR  + UNR 

where;  UTR  = (UR*v) v 

and,  UNR  = UR  - UTR 

From  UNR  and  UTR,  normal  (DN)  and  tangential 
(DT)  drag  forces  can  be  calculated  using  the  pertinent 
drag  coefficients  (C^  and  C^)  . Representative  areas  of 
the  differential  element  ' ds'  can  be  calculated  using  the 
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reduced  diameter  and  stretched  length.  These  are  given  by 
dAj,  and  dA^.  The  standard  formulation  is  given  by: 

DNds  = r/2  (CDt-)  (dA^)  I UITr  | UNR 

and  DTds  = c/2  (CDT)  (dAT)  | UTR  j UTR 

The  constant  force  due  to  gravitational  attraction 

A 

is  W = - W ds  i2.  Here,  is  the  weight  in  air  per  unit 
stretched  length  of  the  mooring  line.  W is  resolved  into 
normal  (WN)  and  tangential  (WT)  components  as: 

A A A 

WT  = (W*v)v  = - t^dscostf^v 

and  WN  = W - WT 

= - W ds(i  - v cosifi.,) 

a 7 *5 

The  hydrostatic  pressure  force  on  the  differential 
element  by  the  surrounding  fluid  is  given  by  the  weight  of 
the  fluid  displaced  minus  the  hydrostatic  pressure  forces 
on  the  end  cross-sections  of  the  element.  Or, 

FHds  = (Wa  - Ww } ds  iz  - *r2[v(Pb  - Pfc)  - 
2 Is  ds  (Pb  + Pt) ^ 


4 


where,  W is  the  weight  in  water  per  unit  stretched  length, 
and  r is  the  reduced  radius  due  to  stretch  of  the  mooring 
line.  and  are  the  hydrostatic  pressures  at  the 

bottom  and  top  ends  of  the  differential  element.  Assume; 


(pb  ~ Pfc)  = eg  ds  cos<J>3 


P + P 
d t. 


2 


P 


c 


and 


W,  - W, 


w 


2 

g frr 


.*.  FH  ds  = (Wa 
- 1 
II. 

iig  3 s 


Ww)ds 


- V COS^3) 


+ 


Here,  Pc  is  the  hydrostatic  pressure  at  the 
midpoint  c of  the  differential  element. 

For  a continuous  cylindrical  body  it  is  also 
assumed  that  the  body  motion  accelerates  water  only  in 
the  direction  normal  to  its  longitudinal  axis.  The 
inertia  forces  due  to  body's  (differential  element) 
own  acceleration  and  the  so  called  "added  mass"  term  can 


now  be  written  as: 


__  wa  - (Wa  - V ~ 

IF  ds  ds  Y — ds  YN 


where;  YN  = Y - YT 

• * ••  A A 

YT  = (Y • v) v;  and  the  added  mass  of  the  differential 
element  is  assumed  to  be  the  mass  of  water  displaced  by 
this  element. 

Combining  ail  these  forces,  the  force  equilibrium 
is  written  as: 


3T  ds 

3v  ds 

3T  ds 

3 v ds' 

T + 

v + — — 

- 

T — — — 

v — * - — 

3s  2 

3s  2 ( 

3s  2 j 

3 s 2 j 

DM  ds  + DT  ds  - Wa  ds  COS^V  - 

A A 

Wa  ds(iz  - v cose p^)  + 

A A 

(Wa  - Ww)ds{  (i2  - V cos^t ) + 


(2.7) 


Or, 


P 3v 

c _ 

eg  3s 


wd  _n.  (w  - Ww) 

ds  (YN  + YT) = — ds  YN  = 0 


3T 

3s 


W. 


= (-DT  + Wfl  costf^  + — YT) v 

g 


(2.8) 


J Pc(Wa-Ww)|3v  _ 

jT  + > — = - DN  + Ww(i  - v coSff  2 ) + 

9 


|2Wa  - Ww] 


Also  we  can  write: 


sy  „ 

= v 

3s 


I’are  T,  the  tension,  v,  the  unit  vector  along 
the  mooring  line,  and  Y , the  geometric  displacement  vector 
are  the  dependent  variables  of  interest.  Equations  (2.8) 
and  (2.9)  can  be  simplified  if  we  define  a new  variable 
called  the  effective  tension  (Te)  similar  to  the  one 
described  in  GOODMAN  and  BRESLIN  (1976) . 


Te  = T + 


Pc(Wa  - W„l 


Differentiating  Te  with  respect  to  s 


3Te  3T  W.  - W.,  3P 

^ u V/  c 
3s  3s  (.g  3s 


- <wa  ' Vcos4>3 


/ 
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Substituting  (2.11)  and  (2.12)  in  (2.8)  and  (2.9)  we 
obtain: 

3Te  W . . 

= - DT  + Vi  COS^  + — YT 

9_  " J 

s g 

3v  - 

and  Te  — = - DN  + Ww(i2  “ v cosi^) 

J s 

• • 

YN 

Equations  (2.13)  and  (2.14)  are  derived  for  a line 
of  an  arbitrary  stretched  length  and  of  consequent  reduced 
diameter.  Also  we  assume,  as  reasoned  in  GOODMAN*  and 
BRESLIN  (1976)  for  materials  obeying  Hooke's  law, that  the 
effective  tension  and  not  the  actual  tension  controls  the 
extensibility  of  the  mooring  lino.  The  extension  of  the 
mooring  line  is  caused  by  (1)  pulling  on  the  line  due  to 
tension  and  (2)  squeezing  of  the  line  due  to  hydrostatic 
pressure.  Volume  I (CHHABRA,  1973)  plotted  non-linear 
curves  between  tension  and  elongation  of  various  mooring 
cables  and  ropes.  Using  the  effective  tension  instead 
of  the  actual  tension,  mooring  line  stretch  and  the 
consequent  reduced  diameter  from  some  initial  reference 
state  are  found  as  discussed  in  Volume  I. 


(2.13) 


(2.14) 
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The  analysis  also  allows  for  an  arbitrary 
number  of  intermediate  bodies  such  as  sensor  packages 
and  subsurface  floats  inserted  in  the  mooring  line.  For 
this  case  equations  (2.13)  and  (2.14)  are  replaced  by: 


Te 


Te 


n+1 


+ F 


n 


(2.15) 


Here,  F is  the  sumnatior  of  gravitational 
attraction,  fluid,  and  inertia  forces  of  the  body  inserted 
between  n,  the  differential  element  above  the  body,  and 
n+1,  below  it. 

For  cylindrical  packages,  the  viscous  drag 
representation  remains  similar  to  the  one  for  cylindrical 
mooring  line.  For  a spherical  or  any  other  shaped  package 
this  viscous  drag  force  is  computed  as: 


DA  = p/2  CDb  ADb  | UR | UR 

where;  CQb  and  ADb  are  the  appropriate  drag  coefficient 
and  area. 

The  gravitational  force  and  the  hydrostatic 
pressure  force  are  combined  to  give: 


NB 


W. 


wb 


i 

z 
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where;  W ^ is  the  weight  in  water  of  the  inserted  body. 

Inertia  forces  for  cylindrical  packages  are 
changed  slightly  from  that  of  a continuous  cylindrical 
mooring  line.  Here  the  assumption  of  the  body  motion 
accelerating  water  only  in  the  direction  normal  to  its 
longitudinal  axis  is  dropped.  Instead  two  added  mass 
terns;  one  normal  and  the  other  tangential  to  the 
longitudinal  axis  are  used.  For  cylindrical  packages 
inertia  force  is  gi'v.  n by: 


(m.  Y + 


"nb  ™ 


♦ ”tb  YT) 


and  for  any  other  shape  this  force  is: 


IF.  = - <m. 


"’ab1 


Here;  is  the  mass,  mVK  and  mfcb  are  the  normal  and 
tangential  components  of  added  mass  for  cylindrical  bodies 
and  mak  is  the  added  mass  for  any  other  shaped  body. 
Summation  of  gravitational  attraction,  inertia,  and  fluid 
forces  gives  F.  Integration  of  equations  (2.10),  (2.13), 
(2.14)  and  (2.15)  along  the  spatial  coordinate,  s,  of 
the  mooring  line  gives  the  dynamic  equilibrium  of  the 
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mooring  system  at  discrete  times.  At  any  time  these  four 
equations  can  be  solved  for  positions,  Y (t) , of  the 
mooring  system.  Positions  Y are  calculated  at  short 
enough  time  intervals,  Ats,  so  that  velocities  Y and 
accelerations  Y during  these  time  intervals  can  be 
computed  as  piecewise  constant.  Accelerations  are  neglected 
at  tQ  and  t-^,  velocities  are  neglected  at  tg . 

2.2.2  Lumped  Parameter  Formulation 

In  this  formulation,  the  mooring  line  is  reduced 
to  a discretized  dynamic  system  which  is  solved  by  a lumped 
parameter  approach.  All  forces  acting  on  the  mooring  line 
and  the  inserted  packages  are  transferred  to  a fixed 
number  of  nodes  (say  N) . If  each  node  has  q degrees  of 
freedom  then  N x q simultaneous  differential  equations,  q 
for  each  node,  describing  these  nodes  can  be  solved 
simultaneously  by  using  any  of  the  numerical  integration 
techniques.  Once  again  the  forces  transferred  at  each 
node  (lumped  mass)  are  the  same  as  listed  in  2.2.1,  with 
the  exception  of  Froude-Krylov  forces  which  are  not 
neglected  here.  The  forces  for  the  ntVl  mass  for  a N node 
system  are  derived  next. 

Consider  the  discretized  dynamic  system  as 
shown  in  Figure  2.6,  where  an  N-mass  system  is  shown.  In 

A 

this  figure  vn  is  a unit  vector  along  the  n*--^  segment 
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(below  the  nts-  mass)  . This  can  be  written  as: 

A A /A  «A 

vn  = ix  coshn  + Ay  cos*2n  + cos*3n 

? , TI+^are  the  tension  vectors,  and  Un  is  the  fluid 
velocity  vector  at  the  location  of  the  nt^1  mass.  A 
typical  mass  could  consist  of  cylindrical  components, 
spherical  components,  and  any  other  shaped  components. 

The  transfer  of  all  adjacent  forces  to  a typical  nfc^  mass 
is  done  as  explained  below. 

Viscous  drag  forces  acting  on  the  mooring  system 
have  to  be  transferred  at  'H'  nodes.  Any  scheme  devised 
for  this  transfer  should  take  into  account  the  different 
behavior  of  cylindrical  and  spherical  components  of  the 
mooring  system.  Let  DNn  be  the  drag  force  transferred 
to  nt:'  node  which  is  equivalent  to  contribution  of  normal 
drag  forces  on  cylindrical  components  lumped  at  the  n4^ 
node.  Similarly,  DTn  is  the  drag  force  on  ntV|  node 
equivalent  to  tangential  drag  forces  on  cylindrical 
components  lumped  at  the  n^*1  node.  DAn  is  the  drag  force 
acting  on  the  n*"^1  node  in  the  direction  of  the  relative 
velocity  URn,  and  is  a contribution  from  spherical  or 
any  other  shaped  components  lumped  4t  the  nt,h  node.  These 

I 

forces  are  given  by  : 


' / 
' / 
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DNn  = (CDIN)nlUNRn|  UNRn 

DT  = (CDIT)  iUTR  ! UTR„ 
n n n 1 n 

DA  = (CDIA)  | UR | UR_ 
n n 1 n 1 n 

and  DF  = DiT  + DT  + DA  (2.16) 

n n n n 

Here;  Un  is  computed  from  the  current  profile  and  tne 
velocity  of  water  due  to  surface  waves  at  the  location 
of  the  nth  mass  which  is  exponentially  attenuated 
with  depth.  CDIN,  CDIT,  and  CDIA  are  the  appropriate 
drag  coefficients  multiplied  by  respective  areas. 

The  gravitational  forces  and  the  hydrostatic 
pressure  forces  are  combined  for  all  the  transferred 
components  to  give; 

NBn  = - Wwn  iz  (2.17) 

Inertia  forces  for  the  components, in  general, are  given 
by: 


IF  =••  - (m  + m ) Y„  - m YN 
n n an  n nn  n 


(2.18) 


- m YT 
tn  n 
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This  representation  is  similar  to  the  one  used 

for  inserted  bodies  in  the  continuous  line  formulation. 

The  tensile  forces  acting  on  the  n*"*1  mass  change 

in  direction  as  well  as  magnitude  with  the  motion  of  the 

nodes.  The  direction  can  be  obtained  by  Y„  and  Y , 

n n-i 

(refer  Figure  2.A);  but  as  the  present  problem  is  one 
of  large  displacements  the  magnitude  is  determined  at  each 
integration  step  by  updating  its  value  at  the  beginning 
of  the  time  step  by  the  stiffness  force  due  to  the 
incremental  deformation  during  this  time  step.  Let  the 
increment a 1 change  in  length  of  the  nth  segment  ALn  be: 


If  kn  (stiffness  coefficient)  is  defined  as  the 

A 

force  required  at  node  n for  a unit  displacement  along  vn; 
then  the  change  in  tension  magnitude  ATfi  is: 

AT  = k AL 
n n n 
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or;  Tn  = (Tn  - ATn)vn  + ATnVn 

Also,  as  tension  vectors  Tn  and  Tn+1  are  acting 
at  the  nth  mass,  we  have: 


TFn  = T 


n+1 


- T_ 


n 


(2.19) 


Froude-Krylov  force  on  the  submerged  'mass  n'  is 
assumed  to  be  given  by  the  product  of  mass  of  water 
displaced  by  the  ' mass'  and  the  acceleration  of  water 
at  that  location  which  is  exponentially  attenuated  with 
depth.  Or, 

_ |K|  zn  rc 

FFn  = mdn  ^5oe  [ixKxcos(K-S-wt)  + 

A A 

i K cos(K-r-u)t)  + i | K 1 sin  (K*S-u>t)  (2.20) 

y y z 

where,  mdn  is  the  mass  of  water  displaced  by  the  mass. 

Combining  equations  (2.16)  through  (2.20)  for  all 
nodes  and  taking  into  account  all  the  boundary  conditions, 
the  force  balance  reduces  in  matrix  notation  to: 


[M]{Y}  = { DF } + {CF } + {TF }+[K]{AY}+{FF} 


(2.21) 


whore; 


[m]  is  a matrix  of  appropriate  inertia  terms. 

[k]  is  a matrix  of  stiffness  terms. 

{TF}  is  an  array  of  tension  components  at  the  nodes. 

{DF}  is  an  array  of  appropriate  viscous  drag  terms. 

{CF}  is  an  array  of  constant  nodal  forces 
(gravitational  and  hydrostatic) . 

{FF}  is  an  array  of  Froude-Krylov  forces. 

(Y)  is  an  array  of  accelerations  at  the  nodes. 

{AY)  is  an  array  of  incremental  deformations  of 
the  nodes. 


derivation  of  Matrices:  In  equation  (2.21),  the 

order  of  each  matrix  and  array  is  equal  to  the  number  of 
nodes  (N)  multiplied  by  degrees  of  freedom  (q  = 3)  per 
node.  For  simplicity  the  masses  were  assumed  to  be  lumped 
at  ' N*  nodes,  in  which  case  the  off-diagonal  terms  in  the 
mass  matrix  are  zero,  i.e.,  force  at  node  ' i'  due  to  an 
acceleration  at  node  * j * equals  zero.  A consistent  mass 
matrix  analysis  can  be  done  to  allow  for  distribute.!  mass. 


Mass  Matrix  - [l*J 


[M]  = 
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where;  M^,  M2,  ....  Mjj  are  3x3  matrices  given  by: 


■ 2, 

m +ra  +m  sin  ©, 
n an  nn  Tin 


+m.  cos  d> 
tn  Yln 


(m  -m__)cos<j),  cosdi  I (in  -m  )cosd>,  cos*., 
tn  nn  In  T2n  1 tn  nn  rln  3n 


j 2 | 

M = (m^  -m  )cos<)).  cos <p,  m +m  +m  sin  <p . I (m  -m  )cosd>„  cos*, 

n tn  nn  2n  In  , n an  nn  T2n  I tn  nn  Y2n 

1 2,  t 

, +mtnCOS  ^2n  | 

f 

I • 2 

(n’tn-n,nn,C0^3nC°^ln!  (mtn'mnn)COs4,3nCOS,*,2n  ; VW*"  *3n 

! ' 2 A 

! , +mtncos  ^3n 


and  m ,m. 

nn  tn 


mass  lumped  at  the  nth  node. 

added  mass  at  n^*1  node  due  to  non- 
cylindrical  components 

Normal  and  tangential  components  of 
added  mass  at  n-h  node  due  to 
cylindrical  components 


Stiffness  Matrix  [k]  - 4 Mass  System 

lr*i . .2  t 2 1 


[K]  = \ 1 


K2  + K3  | 


K3  + K4 


0*0  I 


where;  K 


and 


matrices 


(TF)  = { 
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are  3x3  matrices  whose  elements  are  given  by 


k cs . cs 
s 1 , s 


j/S 


ks  = Stiffness  coefficient  for  segment  s. 

cs . = cosi>.  of  segment  s. 

X / S X 

{TF},  {DF},  { CF } , {FF},  {Y}  and  {AY;  are  column 
These  matrices  can  be  written  as 


(Tj  - AT1)  CEliX  + (T;  - tt T2)  CSx  , 

FFl,l' 

<T1  - AT1>  CS2,1  + (T2  - AT2>  CS2,2 

FF2 , 1 

<T1  - AT1>  CS3,1  + (T2  - AT2>  CS3,2 

FF3 , 1 

• 

> { FF } = J 

1 

(TN-1  “ ATN-1)CE3,H-1 + (Tn ~ ATN>CE3,N 

. 

- <T»  - ATn)csi,n 

FF 

1 , N 

- (T„  - atn)cs2iN 

FF2,N 

- «„  - aVcs3,n 

FF3 , N 

X 
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DF1,1 

CFl,l 

Yl,l 

AY1,1 

DF2,1 

CF 

<■  1 1 

Y2,l 

AY2,1 

DF3,1 

CF3 , 1 

Y3,l 

AY3, 1 

• 

• 

■;  {cf}  = • 

• 

• 

. 

•;  {Y}  = ■ 

• 

• 

•?  {AY}  - • 

• 

. 

• 

dfi,n 

• 

CF1,H 

• 

• • 

Y1,N 

• 

AY1,N 

DF2,N 

CF2,N 

•• 

V 

2,N 

AY2,N 

DF3,N 

CF3,N 

- 

4 

AY3 , N 

Equation  (2.21)  can  be  re-written  as 

(Y)  = [K]"1  {TF}  + [M]'1  [Kj  {AY}  + 

[M]"1  {DF}  + [M]_1  {FF}  + [M]"1  {CF} 


These  3N  second-order  differential  equations  can  be 
reduced  to  a set  of  6N  first-order  differential  equations 
which  can  be  solved  by  any  of  the  numerical  methods. 


'v 

\ 
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2 . 3 Attachment  Between  Surface  Buoy  and  a Mooring 

Line 

Figure  2.7  shows  the  attachment  between  a surface 
buoy  and  a mooring  line.  Let  the  attachment  point  be  'a', 
which  is  fixed  rigidly  to  the  buoy. 


-T 

Figure  2.7  Attachment  Between  Surface 
Buoy  and  a Mooring  Line 

In  addition  to  all  the  other  forces  acting  on  the 
surface  buoy  as  outlined  in  Section  2.1;  a tension  vector 
T is  acting  at  point  a,  which  should  be  included  in  the  force 
balance . 

T = Tv  = Txix  + Tyiy  + Tziz  (2.22) 


/ 
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also  let  vector  Rca  be  given  by: 


Rca  * lx'  xa*  + V Ya'  + V V 


Taking  moment  of  T about  point  c we  obtain: 


TM  = - R x T 
ca 


for  inclusion  in  the  moment  balance. 


Two  Dimensional  Analysis; 


Rca  - V xa'  + V *a' 


= i„(x='cos0  + z 'sin6)  + i„(z  'cos0  - x_ 

X a a Z a a 


and; 


TM  = i. ,[T„  (x=  'cos6  + z'sinQ)  - T (z  'cos0 

y Z a a X a 


Three  Dimensional  Analysis: 


Rca  = Ax  (xa  ” ^a  + ®za)  + V*xa  + *a 


♦*;) 


iz(-ex;  + + *;> 


(2.23) 

' sin0 ) 
x, ' sir.0)  ] 

O 


I -V «»/*.. 
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And;  TM  = ixLTy  (-8x^  + 4>y  'p  + z'a)  - Tz  - <|)Z^)  ] + 

iy[Tz(x4_^a  + 0za)  ' Tx  (_9x4  + + 2a}  ] + 

iz[Tx(^xa  + yk~  *za}  ‘ Ty(xa‘  ^k+  6z4^ 


Change  in  mooring  line  tension  magnitude  due  to  buoy  rotation 
(coefficients  in  the  stiffness  matrix) ♦ 


Two  Dimensional  Analysis: 


Rotate  R 

ca 

through 

A0  to 

obtain 

R«i 

1 

A9 

cose 

sin0~J 

kl 

V 

L-40 

1 

|j-sine 

1 

cosej 

< 

zk 

For  small  rotations  change  in  tension  magnitude  AT  is  given 
by 


AT 


= k (R 


ca^,  Rca^ 


. v 


where  k is  the  stiffness  coefficient  for  the  line  segment 
attached  to  the  surface  buoy. 

A A A 

If  v = ixcos<{>  + iz  sin$ 


AT  = k[cos  t'  (Zg  cosO  - x^  sinO)  + simM-z^  sin0  - 
x'  cos9)]A6 

a 
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Three  Dimensional  Analysis; 


K 


ca 


1 


r - 

| 1 

1 o 

'‘a 

v'  1 -A* 

<l>  l 

y'a 

A>  1_ 

-!)  f 1 

zl  ! 
aJ 

and; 


AT=k 


cos  :2(  Ox;-  :y;-z;)!cos ; x (-.•x;+^y;+z;);cost  x (-'ix;-y;+i  z;) 

L.  L._ 

I I 

+005*3  (+x;+y;-'{iz;)j+cos?  3 (-x;+>j/y;-oz;)!+cos4>2  (x;-^y;+ez;) 


i * . 
i • . . 

i 


If  more  than  one  mooring  lines  are  attached  to  the 
surface  buoy,  as  is  the  case  of  simulation  study  in  Section  5.3, 
force  and  moment  of  both  these  lines  are  to  be  included 
in  the  appropriate  equations.  It  is  important  to  note  that 
within  the  definition  of  the  coordinate  frame  (fixed  at 
the  anchor)  if  a mooring  line  hangs  from  the  buoy  (as  is  the 
case  in  5.3),  - T of  this  section  would  be  replaced  by  + T. 

Also  in  this  case  AT  would  be  given  by 


= - k ( R 


ca. 


R ) • v 
ca 


AT 


2 . 4 Window  Shade  Drogue 

Sometimes  an  alternative  approach  to  mooring 
surface  buoys  is  the  use  of  surface  trackable  drogued 
drifting  buoys.  Such  a buoy  system  employs  a high  drag 
device  (or  drogue)  tethered  to  a trackable  buoy  at  the 
surface.  A drogue  is  subjected  to  the  sane  forces 
(including  the  wave  effects)  as  a differential  element 
of  tne  mooring  line.  A special  case  of  drogues,  the 
window  shade  drogue,  as  shown  in  Figure  2.3  ie  considered 
in  this  analysis.  The  drogue,  modeled  dynamically  as 
a lumped  parameter  system  (Figure  2.8),  is  reduced  to 
strips.  The  mass,  the  acting  forces,  and  the  dynamic 
behavior  of  each  strip  is  lumped  at  an  assigned  nodal 
point.  Nodal  masses  are  linked  by  suitably  elastic  lines 
The  window  shade  drogue  is  a compliant  sheet. 
Close  observation  of  scale  model  drogues  (Vachon,  1973) 
has  revealed  that  it  seems  nearly  locked  to  the  local 
water  mass  insofar  as  motion  normal  to  the  '•heel,  is 
concerned.  However,  in  the  tangential  direction, it  seems 
to  slip  through  the  local  water  mass  rathei  easily.  The 
forces  were  therefore  described  as  follows:  the  wave 

forces  (Froude-Kry lov  force;  as  it  is  assumed  that  there 
is  no  diffracted  wave)  acting  on  a material  strip  are 
estimated  by  calculating  the  force  normal  to  the  axis 


Surface  Buoy 


Figure  2.8  Drifting  Drogued  Buoy  System 
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of  revolution  of  a hypothetical  cylindrical  element,  the 
diameter  and  length  of  which  is  equal  to  the  breadth  and 
length  of  the  strip  respectively.  This  force  is  given 
by  the  mass  of  the  hypothetical  water  cylinder  multiplied 
by  the  component  of  water  acceleration  normal  to  the 
cylinder' s longitudinal  axis.  Or; 

""dn'^Ew  - <SEW''n>''n5  <2-24 

where  w is  the  particle  of  water  next  to  the  n*"*1  node, 
and  m^n  is  the  mass  of  the  hypothetical  water  cylinder. 

For  inertia  forces  man  and  mfcn  equal  zero;  whereas  mnn 
equals  ndn.  The  viscous  drag  of  this  strip  in  the  normal 
direction  is  calculated  using  the  frontal  area  of  the  strip 
and  the  drag  coefficient  measured  in  water- filled  quarry 
tests  (Vachon,  1975) . The  elemental  strips  of  the  drogue 
are  assumed  subject  to  tangential  viscous  drag  calculated 
using  the  area  of  the  strip  and  a tangential  drag 
coefficient  measured  during  drop  tests  in  the  ocean 
(Vachon,  1975) . 

Gravitational,  hydrostatic  pressure,  and  tensile 
forces  on  a drogue  strip  are  similar  to  the  ones  outlined 
in  Section  2.2  for  a discretized  dynamic  system. 
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3.0  METHOD  OF  SOLUTION 

Section  2 outlined  all  the  necessary  ingredients 
of  the  mathematical  models  used  in  this  report.  In  this 
section  some  details  of  how  to  use  these  ingredients,  to 
obtain  a solution  for  a specific  problem  on  hand,  are 
given.  For  illustrative  purposes  a specific  case  study, 
of  a spar  buoy  at  the  surface  and  a mooring  line  connecting 
this  buoy  with  the  anchor,  will  be  discussed.  Method 
of  solution  for  both  the  two-dimensional  and  the  thrae- 
dimensional  analysis  of  the  moored  system  will  be  outlined  in 
3.1.  Section  3.2  outlines  the  approximate  initial 
conditions  for  the  steady  state  analysis  of  a spar  buoy 
freely  floating  in  a surface  wave. 

Another  case  study  of  a spar  buoy  at  the  surface 
and  connected  to  a window  shade  drogue  at  a given  depth 
is  discussed  in  section  3.3. 

3. 1 Moored  System  Analysis 

Two-Dimensional; 

In  this  case,  the  current  profile  is  in  the  same 
plane  as  the  train  of  surface  gravity  waves  having  a 
single  frequency  and  amplitude;  ana  its  direction  parallel 
or  anti-parallel  to  the  direction  of  surface  wave 
propagation. 
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Static  Solution: 

For  the  static  solution,  the  analysis  is  started 
at  the  spar  buoy.  The  pertinent  static  equations  for  the 
soar  buoy  (Section  2.1)  with  an  attached  mooring  line 
(Section  2.3)  are: 

DFx  " Tx  = 0 (3.1) 

Pog  - Mg  + DFz  - Tz  = 0 (3.2) 

and  -Pig  sinO  + DM  + T_  (x'  cosO  + z*  sin0)  - 

X Z a a 

- T ( z ' cosO  - x‘  sin0)  = 0 (3.3) 

X a a 

Here;  we  have  three  equations  and  four 
(T  , T , 6,  h ) unknowns,  where  h is  the  static  draft  of 

X Z 9 b 

the  spar  buoy.  Now  if  we  assume  a value  for  hs,  the 
other  three  unknowns  can  be  solved  by  equations (3. 1) 
through  (3.3).  The  solution  for  the  three  unknowns  T , 

Tz,  and  P is  found  by  iterating  on  the  pitch  angle  0. 

If  0 is  known;  DFX,  DFz,  and  Dfl  can  be  computed  explicitly. 
PQ  and  P^  can  be  calculated  if  hg  is  known  and  thus  Tx 
and  Tz  are  given  by  equations  (3.1)  and  (3.2)  respectively. 
Moment  balance  can  now  be  performed  by  using  equation  (3.3) 


I 
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to  obtain  the  error.  This  ei ror  is  minimized  by  iterating 
on  the  pitch  angle  0. 

Having  solved  the  buoy  force  and  moment  balances 
for  a assumed  hg  the  mooring  line  can  be  analyzed  next. 

For  the  continuous  line  formulation,  a numerical 
integration  is  carried  down  the  mooring  line  by  dividing 
it  into  a finite  number  of  segments.  In  this  representation 
equations  (2.13),  (2.14),  and  (2.10)  for  two-dimensions 
can  be  written  as: 


ATe  = (-DT  + Ww  sinct>)  AS 


(3.4) 


TeA<*>  = [DN  + W cos<t>]AS 


(3.5) 


AY  = i AS  cosiji  + iz  AS  sin<f> 


(3.6) 


Having  found  T„  and  Tz  (or  T and  0)  at  the  top 
of  the  mooring  line  equations  (3.4)  through  (3.6)  are 
used  to  find  the  configuration  and  the  tension  distribution 
along  the  mooring  line.  At  the  end  of  the  mooring  line 
the  boundary  condition  of  the  ocean  floor  depth  should 
be  net.  If  not,  the  whole  process  is  repeated  with  a new 
h_;  thus  computing  the  correct  value  of  h_  for  the  given 
bound.' ry  condition.  Whenever  the  mooring  line  is 
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discontinued  to  insort  an  intermediate  body,  the 
discontinuity  is  taken  as  a lumped  mass  and  resolved  as 
done  below  for  the  lumped  parameter  approach. 

Tor  the  lumped  parameter  f ormulation,  the 
configuration  ar/’  the  tension  distribution  along  the 
mooring  lino  can  be  found  by  (Section  2.2): 


= T 


n+1 


W 


wn 


+ DF 


n 


(3.7) 


where;  T is  the  tension  vectot  below  the  n1  mass  and 

1 , . abcve  it. 
n+1 

Dynamic  Solution: 

Equations  of  Section  2.2  which  are  written  for 
three-dimensional  analysis  can  be  easily  reduced  for  the 
two  dimensional  problem.  Using  the  lumped  parameter 
approach  for  the  mooring  line;  equations  (2.21)  are 
combined  with  equations  (2.3),  (2.22),  and  (2.23)  to 

obtain  a global  matrix  equation  for  the  mooring/buoy 
system.  Stiffness  coefficients  for  the  change  in  mooring 
line  tension  magnitude  due  to  buoy  rotation  are  obtained 
from  Section  2.3.  This  matrix  equation  can  be  solved 
in  the  time  domain  using  various  numerical  integration 
techniques.  Integrals  of  and  given  in  Section  2.1, 


/ 


\ 


\ 


/ 


r * 
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are  computed  exactly;  whereas  the  integrals  for  viscous 
drag  forces  and  moments  DFX,  DFZ,  and  DM  are  computed 
numerically  by  dividing  the  spar  buoy  in  disks  of  height 
dz'.  The  composit  matrices  are  of  order  2N  + 3.  Matrices 
[M],  (TF),  (K),  (DF),  (FF),  {FD},  {FG},  and  (FA)  are 
computed  at  each  time  step  of  the  integration. 

For  the  continuous  line  formulation;  the  matrices 

of  equation  (2.3)  combined  with  equations  (2.22)  and  (2.23) 

can  be  solved  for  T in  equation  (2.22)  if  the  buoy  draft 

hQ  is  known.  Then  integration  of  equations  (2.13), 

(2.14),  (2.10),  and  (2.15)  along  the  spatial  coordinate,  s, 

of  the  mooring  line  gives  the  dynamic  equilibrium  of  the 

mooring  system  at  discrete  times.  Iteration  on  the 

unknown  h may  be  necessary  to  meet  the  boundary  condition 
o 

of  the  ocean  floor  depth. 

Three  Dimensional; 

Here;  the  train  of  surface  gravity  waves  having 
a single  frequency  and  amplitude,  can  propagate  in  a 
general  direction.  The  current  profile  is  also  three- 
dimensional  and  can  vary,  in  magnitude  and  direction,  with 
depth. 
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Static  Solution; 

The  pertinent  static  equations  for  the  spar  buoy 
(Section  2.1)  with  an  attached  mooring  line  (Section  2.3) 
are: 


and; 


) 


dB(Rcp  x g) 


+ DF  - T = 0 


+ DM  + TM  = 


0 


(3.8) 

(3.9) 


Integrations  for  dB  are  over  the  static  draft  of  the 
buoy,  and  integration  for  DM  is  over  the  length  of  the  buoy. 
Equations  (3.8)  and  (3.9)  can  be  written  as: 


Pog  iz+DF-T=0 


(3.10) 


and 


e 


Ls]- 1 (F) 


l* 


(3.11) 


where;  [s]  is  a 3 x 3 matrix  with  components  given  by: 


’ii  - VH  + dV2  - Vi  ' Vi 
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and 


also 


12 


Ty'xa  » 


’13 


T x' 

z a 


’21 


T y ' ; 
xJa  ' 


23 


zJa 


S22  * V'l  + °?zl2  - Vi  - Vi 


S31  * - DFxl2  + Vi 


S32  - - DFyl2  + V* 


S33  - ' Vi  ' Tyya 


{F}  is  a 3 x 1 array 


F1  ‘ - DFy*2  + Tyza  - Vi 


F,  = 


DFxJ2  - Vi  + Vi 


F3 


Vi  - Tyxi 


dB;  and  £.2  are  the  distances  between 


center  of  gravity  and  center  of  buoyancy  and  center  of 


drag  respectively.  For  the  static  case  we  also  assume 
(small  angles) : 


DF  « V V 
x 'Bwx1  Bwx 


DF  « I V V 
u y ‘Bwy1  Bwy 


and 


DF 


V V 
1 Bwz1  Bwz 


.*.  T and  <$> , 0,  and  ip  can  be  solved  explicitly  from 
equations  (3.10)  and  (3.11)  if  hg  is  known.  hg  is 
determined  by  iteration  on  the  boundary  condition  of  the 
ocean  floor  depth. 

Mooring  line  solution  is  similar  to  the  one 
presented  for  the  two-dimensional  analysis. 

Dynamic  Solution; 

This  dynamic  solution  is  similar  to  the  one  fo 
two-dimensional  analysis  except  that  the  order  of 
matrices  in  this  case  is  3N  + 6. 


3. 2 Initial  Conditions  for  the  Steady-State 

Analysis  of  a Spar  Buoy 

Two-Dimensional  Analysis: 

Let  CN  = 1.0,  CT  = 0.0,  and  xc  = 0.0.  Also 
assuming  small  pitch  of  the  buoy  and  neglecting  viscous 
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drag  forces,  the  equations  of  motion  for  a spar  buoy 
become : 

2Mxc  + P^'s  = 2<j 2£0Qq  coswt  (3.12) 

••  2 
M*z„  + oS  _gz  = - w f Q„  sinwt 
c o2  o o 

+ pSo2gCQ  sinwt  (3.13) 

Pxxc  + (Iyy  + P2)0  + P1g6  = 2w2^oQ1  coswt  (3.14) 
M*  **  M + 4/3  p[r^3  + (rj-r-,)3] 

Surge  Equation  3.12: 

Let  P^'o  <<  2Mxc  and  x = C^sinwt  + C2  coswt  + C^. 
Now  differentiating  x twice  and  substituting  in  (3.12)  we 
obtain 

x = - ?0Q0  coswt/M  + C3 
and  x = w^0Q0  sinwt/M  + C^t 

Hence  x(o)  = 0 


and 


where , 


= - ( Q /H 
'o  o 


and 


x (o) 
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Heave  equation  (3.13) 

Let  7 - sinwt  + C2  coswt  + C3 

Differentiating  twice  and  substituting  in  equation  (3.13)  we 
obtain: 


C2  = c3  = ° 

and  = C0(pSo2g  - oj2Qq)  / (pSQ2g  - M*w2) 

Hence,  z (o)  =0 

z(o)  = ujC^ 

Pitch  Equation  3.14: 

Substituting  xc  from  equation  (3.12)  into 
equation  (3.14)  we  obtain:  « 

A0  + B0  = C coscot  (3.15) 

where,  A = (I  + P2  - PX2/2M) 


B = Pxg 
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and  C = (2<d2/,0Q1  - PjW^Qq/M) 

Now  let  B = sinwt  + coswt  + 

Differentiating  0 twice  and  substituting  in  equation 
(3.15)  gives: 

C.  = C-,  = 0 
1 

and  C0  = C/(B-Aw2) 

c 

Or,  0 = s-  coswt 

B-Aw 

Hence, 

. " «o(2Oi-piVM) 

0 (o)  — ' nj  ~ ^ 

Plg~u'"(1yy+PZ'Pl  /2M) 

and  0 (o)  = 0 

/ 

Three  Dimensional  Analysis: 

Similarly;  for  the  three  dimensional  case  with 


wave  direction  at  an  angle  0 with  the  x-axis  we  have: 
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x(o)  = - Kq  cosB  Qq/M 
y(°)  = - K0  sinB  Q0/M 
z(o)  = 0. 
x(o)  = y (o)  = 0 


% u5o(pSo29~“2qo) 

2(0)  = 2 

pSo2g-M*w 

*} 

w^C.sinB (2Q.-P.Q  /M) 

•Mo)  2-= ±—±  ° 

P1g-u.z(Iyy+P2-P1V2M) 

w2C0  cosB (2Q1-P1Qo/M) 

e (0)  = 

pi dyy+P2-P1V2M) 


ip  (o)  = 0 


4>  (o)  =0  (o)  = ip  (o)  = 0. 


3. 3 Drifting  Drogued  Spar  Buoy 

\ 

In  this  analysis  the  drogue  and  the  mooring  line 
connecting  it  to  the  spar  buoy  are  modeled  dynamically 
as  a lumped  parameter  system. 
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Two-Dimensional  Analysis 
Static  Solution: 

The  pertinent  static  equations  for  the  spar  buoy 
(Section  2.1)  with  an  attached  mooring  line  (Section  2.3) 
are: 


DFx  + Tx  = 0 (3.16) 

Pa  - Mg  + DF  + T_  = 0 (3.17) 

oJ  J z z 

and  -P.g  sin0  + DM  - T_  (x'  cos8  + z!  sin0) 

x z <1  a 

+ T (z'  cos0  - x'  sin0)  = 0 (3.18) 

x a a 

Here;  we  have  three  equations  and  five  (V  , h . 

b S 

T , T , 0)  unknowns  where  V is  the  static  velocity  of 
the  drogued  drifting  system.  These  equations  are  solved 
by  iterating  on  V and  0.  Assuming  a value  for  V , T 
and  Tz  can  be  calculated  by  known  weights,  buoyancies, 
viscous  drag  forces,  and  elastic  characteristics  of  the 
system  beneath  the  spar.  Now  if  a value  for  0 is  assumed 
equation  (3.17)  is  solved  for  hs  and  the  error  in 
equation  (3.18)  is  found.  This  error  is  minimized  by 
iterations  on  9.  Next  Vg  is  iterated  upon  to  satisfy 
equation  (3J.6)  . 
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Dynamic  Solution 

This  solution  is  similar  to  the  lumped  parameter 
solution  discussed  in  Section  3.1. 

Three  Dimensional  Analysis 
Static  Solution: 

Here  the  pertinent  equations  (static)  remain  the 
same  as  Section  3.1  expect  for  replacing  -T  by  +T  and  + TM 
by  -TM. 

Also;  following  similar  reasoning  of  small  angles 

T,  <)>,  0,  and  ip  can  be  solved  explicitly  from  equations  (3.10) 

and  (3.11)  with  new  T and  TM  if  V„  is  known.  V is 

s s 

determined  by  iteration  on  the  equation  similar  to  (3.10). 

Dynamic  Solution 

This  solution  is  similar  to  the  one  in  Section  3.1  . 

\ 

\ 

\ 


/ 


/ 
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4.0  COMPUTER  PROGRAM  DETAILS 

This  section  will  describe  the  various  computer 
programs  written  for  the  theoretical  analysis  presented  in 
Section  2 and  using  solution  methods  of  Section  3. 

Description  of  the  programs  is  general  and  specific  details 
can  be  found  in  the  program  listings  presented  in  Appendix  A. 

All  computer  programs  are  written  in  FORTRAN  IV  and  have 
been  run  for  many  case  studies  on  the  CSDL  AMDAHL  470  V6 
computer.  Inout  data  required  for  these  programs  is  also 
explained  in  this  section.  Some  of  the  simulations,  along 
with  the  accompanying  input  data,  are  presented  in  Section  5. 

Computer  programs  written  for  tne  surface  moored 
and  drifting  systems  are  tabulated  in  Table  4.1.  Programs 
for  subsurface  moored  systems  are  tabulated  in  Table  4.2. 

Program  listings  for  some  of  these  computer  programs  are 
presented  in  Appendix  A.  Listings  not  preserted  are  duplications 
of  what  is  presented  and  can  be  obtained  after  minor  modifications. 

4.1  Surface  Moored/Drifting  Systems 

Computer  program  SD3.FORT  is  coded  for  the 
three-dimensional  analysis  of  surface  moored  and  drifting 
systems  presented  in  Section  2.  Four  case  studies  can  be 
simulated  by  this  program: 


1.  A freely  floating  spar  buoy. 


FILE  NAME  MOORING  SYSTEM  DEGREES  OF  FREEDOM  PHYSICAL  FORCES  j INPUT  ENVIRONMENT 

I (Forcina  Function) 
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TABLE  4.1  Programs  for  Surface  Moored/Drifting  Systems 


TABLE  4 
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2.  A spar  buoy  anchored  with  a lumped 
parameter  mooring  line. 

3.  A spar  buoy  anchored  with  a lumped 
parameter  mooring  line  and  a lumped 
parameter  mooring  line  hanging  from 
the  buoy. 

4.  A spar  buoy  attached  to  a window  shade 
drogue . 

SD2.F0RT  does  the  same  simulations  as  SD3.F0RT, 
but  is  coded  for  the  two  dimensional  analysis  of  Section  2. 
Computer  programs  BUOY. FORT,  MDE.FORT,  and  DDB.FORT  are 
subsets  of  the  computer  program  SD3.F0RT.  BUOY. FORT 
simulates  case  study  1,  MDE.FORT  simulates  case  study  3, 
and  DDB.FORT  simulates  case  study  4.  All  these  programs 
use  the  lumped  parameter  formulation  for  the  mooring  line. 
Input  data  necessary  for  these  programs  is  outlined  below. 

Input  Data  Required  for  Computer  Program 

SD3. FORT  (All  units  in  F.P.S.) 

General  data: 

NM:  Number  of  nodes  (masses) . Equals 

one  for  case  study  1. 

NB:  Number  of  the  surface  buoy.  Starting 

from  the  anchor,  the  number  of  the  node 
(anchor  is  not  counted  as  a node)  where 
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ND: 

NST: 


NC: 

MOOR: 


DEPTH: 

DT: 

TMAX: 

T2: 

Buoy  Data 
RD1, RD2 : 


the  buoy  is  located.  NM  = NB  for 
case  study  3.  NB  equals  one  for 
case  studies  1 and  4. 

Number  of  the  first  node  on  the  drogue. 
Number  of  steps.  This  is  used  for 
dividing  the  spar  buoy  in  NST  strips 
for  integration  of  viscous  drag 
forces  and  moments  over  the  entire 
submerged  buoy. 

Number  of  cycles,  i.e.  number  of 
surface  waves  for  which  the  simulations 
are  tc  be  done. 

Index  for  case  study  control.  Equals 
zero  for  case  study  1;  one  for  case 
studies  2 and  3,  and  two  for  case 
study  4. 

Ocean  depth. 

Integration  time  step  AT. 

Maximum  time  of  simulation. 

Time  interval  for  simulation  printout. 

Two  radii  (r1  and  r2)  of  the  tuned 
spar  buoy.  RDl  is  for  the  drum  and 
equals  RD2  for  a cylindrical  spar. 
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ZCG : 

RGYR: 

HMAX: 

CDL: 

HST: 


CON: 
COT: 
CDP : 
ALPHA: 


Distance  (Z  >0)  between  centeroid  of 
c 

the  spar  base  and  its  C.G. 

Radius  of  gyration  of  the  buoy  about 
any  axis  in  the  horizontal  plane. 
Length  of  the  spar  buoy  (H) . 

Normal  viscous  drag  coefficient  (CDN) 
used  for  the  spar.  ^C^  is  found  by 
multiplying  CDLxO.02 

Length  of  the  drum  for  the  tuned  spar 
buoy  (h^)  . For  cylindrical  spar  IIST 
is  some  arbitrary  number  between  zero 
and  the  submerged  depth  of  the  spar. 
Added  mass  coefficient  (Cjj)  for  the 
buoy.  (Used  & read  in  SD2.FORT  only). 
Added  mass  coefficient  (CT)  for  the 
buoy. (Used  & read  in  SD2.FORT  only). 
Viscous  drag  coefficient  for  the  base 
and  step  of  the  tuned  spar  (CDp) . 
Added  mass  coefficient  for  the  base 
and  step  of  the  tuned  spar  (a). 


Mooring  Line  Data: 

W(I) : Weight  in  water  (w^.)  of  components 

lumped  at  the  I*"*1  node.  Equals  weight 
in  air  for  the  buoy  (NB^  node) 


/ 
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CM(I):  Mass  of  components  lumped  at  the  I ^ 

node  (nr ) . 

CMS (I);  Added  mass  of  spherical  components 
lumped  at  the  node  (m  ■ ) . 

CMD(I):  Mass  of  water  displaced  by  components 

lumped  at  the  node  (n,;). 

C*  X 

CMN'  ( I ) # domal  and  tangential  components  of 
CMT ( I ) : 

added  mass  of  cylindrical  components 

lumped  at  the  I*"*1  node  (m  . and  m.  ■). 
r m ti 

SL(I):  Slack  length  of  the  mooring  line 

proceeding  the  Ith  node  (L^ ) . 

EK ( I ) : Elastic  coefficient  for  SL(I)  ( k ^ ) . 

CDIN(I),  normal  and  tangential 
CDIT(I) : 

viscous  drag  constants  (1/2  pCpAp) 
of  cylindrical  components  lumped  at 
the  I*"*1  node. 

CDIA(I):  Viscous  drag  constant  of  non- 

cylindrical  components  lumped  at  the 
I node . 

DEP(I):  Depth  (less  than  zero)  of  the  I*"*1 

node  (approximate)  (z^)  . 
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Attachment  Data: 

XC1,YC1,  Coordinates  of  the  attachment  point 

r.Cl: 

between  spar  and  the  anchoring  line 
in  spar  buoy  coordinate  frame  (x* , 

y' , z' ) . 

XC2,YC2,  Coordinates  of  the  attachment  point 
ZC2: 

between  the  spar  and  the  instrument 
line  in  x',  y',  z’  frame. 

Current  Profile  Data: 

V(IJ,I):  Absolute  velocity  of  water  in  J 

(x,  y,  and  z)  direction  at  the  depth 
of  node. 

Surface  Wave  Data: 

WE : Frequency  (u)  of  the  surface  wave. 

AMP:  Amplitude  (CQ)  of  the  surface  wave. 

BETA:  Wave  direction  (8). 

All  of  the  input  data  are  read  in  the  main 
program.  The  main  program  calls  subroutines:  (a)  STATIC; 

which  calculates  the  mean  configuration  or  steady  state 
(static)  solution  without  any  surface  wave  forcing. 

(b)  BUOYS;  which  calculates  the  mass  matrix  [ll]  for  the 
spar  buoy.  (c)  DBDRG,  which  calculates  the  viscous  drag 
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forces/monicnts  or!  tha  spar  buoy.  (d)  MINVER,  which  inverts 
a matrix.  (e)  MATR  X,  which  calculates  the  wave  exciting 
plus  hydrostatic  forces  on  the  spar  buoy,  the  stiffness 
matrix  |.K],  the  remaining  mass  matrix  for  the  mooring 
lines,  viscous  draq  forces  on  the  mooring  lines,  Froude 
Krylov  forces  on  the  mooring  lines,  and  the  tension  force 
components  at  each  node.  It  also  in  egrates  the 
differential  equations  of  notion.  (f)  NEXT,  which  updates 
the  geometry  of  the  system. 

Before  finding  the  steady  state  dynamic  response 
of  a system  to  a sinusoidal  surface  wave,  an  equivalent 
current  profile  has  to  be  determined  to  compute  the  mean 
configuration  of  the  system.  This  current  profile  depends 
on  the  original  current  profile  plus  the  contribution  due 
to  the  surface  wave.  In  effect  it  is  a rectification  of 
the  time  constant  current  profile  by  the  oscillatory 
velocity  field  generated  by  the  surface  wave.  From  Section  2.1 
we  have : 1 


'Ew 


= i (V  + <f>  ) 
x ox  x 


+ i (V  + $ ) 
y oy  y 


+ VvoZ  + V 


also; 


I v_  | 2 — (V  + $ )2  + (V  + <i>  ) 2 
1 Ew 1 'ox  x oy  y 


* (Voz  + V 
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And  as  viscous  drag  forces  determine  the  moored 
or  drifting  system  configurations  for  the  static  solution, 
and  these  are  proportional  to  !V£WiV£W»  the  square  roct 
of  mean  ^V£WiV£W)  over  a wave  period  (T)  would  give  the 
equivalent  rectified  current  vector  at  this  location. 


Jr  j |VEwl(VEw»d(“t) 
’ o 


= i E+i  E+i  E 
x x y y z z 


Now  the  equivalent  current  vector  can  be  written  as: 


V = E/EE 
eq 


where:  EE  = ( |E| 


This  calculation  to  find  V.  , before  a static 
solution  is  performed,  is  done  in  the  main  program.  The 
numerical  integration  of  differential  equations  in  the 
subroutine  MATRIX  was  performed  by  the  rectangular  rule, 
which  was  found  adequate  for  all  case  studies  simulated. 


Subsurface  Moored  Systems 


Computer  program  SSD3.F0RT  is  coded  for  the  three 


dimensional  analysis  of  subsurface  moored  systems 
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(continuous  lino  formulation)  given  in  Soction  2.  The 
program  simulates  the  response  of  a subsurface  moored 
system  to  the  time  varying  current  fields.  SSD3.F0RT 
reads  in  the  relative  (relative  to  the  mooring  system) 
velocity  profile  in  the  three  orthogonal  x,  y,  and  z 
directions.  From  these  three  components  the  viscous  drag 
forces  can  be  computed  directly.  These  components  of 
the  relative  velocity  could  be  the  ones  measured  by 
current  sensors  on  the  mooring  line.  Alternatively, 

SSD31. FORT  reads  in  the  absolute  (relative  to  the  -ean 
floor)  velocity  profile  in  the  three  orthogonal  (x,  y, 
and  z)  directions.  Profile  is  read  in  as  amplitudes  of 
the  sinusoidally  varying  currents.  All  components  (varying 
with  depth)  are  driven  by  a single  frequency.  In  addition 
a surface  wave  of  given  amplitude,  frequency,  and  direction 
generates  velocity  field  exponentially  attenuated  in  depth. 
The  program  computes  the  total  absolute  velocity  and  then 
substracts  the  mooring  line  velocity  from  this  absolute 
velocity  to  compute  the  visco-.s  drag  forces.  SSS3.F0RT  is 
the  same  as  SSD3.F0RT  but  it  neglects  all  inertia  forces 
of  the  mooring  system.  All  three  programs  use  the 


continuous  line  formulation  of  the  mooring  line,  given  in 
Section  2,  as  the  mathematical  model.  Innut  data,  except 


for  the  input  forcing  function  (velocity  profiles  and  surface 
wave)  are  the  same  for  all  three  programs  and  aredetailed 
below : 


General 

NIs 


NP: 

NS: 

NPT: 

IKK: 

DEPTH: 

DDT: 

ER: 


ER1 : 


Data: 

Number  of  intermediate  bodies 
(instrument  clusters  or  inserted 
floats)  to  be  included  in  the  analysis. 
Number  of  mooring  parts  (each  mooring 
part  can  have  different  properties). 
Number  of  segments  the  continuous 
line  is  to  be  divided  in. 

Number  of  dynamic  equilibriums  to  be 
performed. 

Number  of  locations  the  velocity  profile 
changes  along  the  ocean  depth. 

Depth  of  the  ocean.  (Meters) 

Approximate  depth  of  the  mooring  line 
top.  (Meters) 

Used  for  iteration  on  the  boundary 
condition  of  ocean  depth.  If  DDT  is 
known  with  reasonable  accuracy,  ER 
is  not  required.  (Meters) 

Used  for  error  control  in  SSD31.F0RT. 


.instrument  in 
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Ir.3*- "UT.ent  Data: 

th 

P I ; • j } Position  of  the  J 

slack  length  distance  from  the 
top.  (Meters) 


SI  (J)  : 


ZM  (J)  : 
ZMC  (J)  : 


ZMV(J)  : 


F (3, J)  : 


CDIN  (J)  , 
CDIT(J)  : 


Length  of  the  Jt^1  instrument. 

Equals  zero  for  non-cylindr ical 
instruments.  (Meters) 

Mass  of  the  J*"*1  instrument  (m^ ) . (Slugs) 

Normal  added  mass  component  (nnj) 

for  cylindrical  instruments.  (Slugs) 

Added  mass  (m  • ) of  the  Jth 
aJ 

instrument  for  a non-cylindrical 
instrument.  Tangential  added  mass 
component  (mtj)  for  cylindrical 
instruments.  (Slugs) 

Weight  in  water  (W.  . ) of  the  J*"*1 
instrument.  (Pounds) 

Normal  and  tangential  viscous 
drag  constants  (1/2  pCpAp)  of  the 
instruments.  CDIT  is  not  used  for 
spherical  instruments.  (F.P.S.  units) 


Mooring  Line  Data: 

th 

DIAL (I) : Nominaldiameter  of  the  a n mooring 


part.  (Inches) 
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SLL(I)  : 


AWL (I)  : 


WWL  ( I ) : 


TPL(I) r 


C01  (I)  , 
C02  (I ) , 
POl(I)  , 
P02(I), 
AO  (I)  : 


CDN(I)  , 
CDT(I)  : 


Current 
D (I ) : 


V(I,J)  : 


til 

Slack  length  of  the  I mooring 
part.  (Meters) 

Weight  in  air  per  unit  slack  length 
of  the  Ifc^  mooring  part  (Wa) . (Jos/M) 

Weight  in  water  per  unit  slacic 
length  of  the  I*"*1  mooring  part  (Ww)  • (los/M) 
Transient  peak  load  on  the  Itn  mooring 
part.  For  materials  which  obey  Kook's 
law  this  is  substituted  with  the  Young's 
modulus  of  elasticity.  (Refer  volume  1 
CHHABRA,  1973).  (lbs.  or  lbs/in2) 

Constants  for  stress-strain 
relationships  of  mooring  line 
materials.  (CHHABRA,  1973), Col  is  used  for 
jacket  diameter  (in.)  in  jacketed  wire  rope. 

Normal  and  tangential  drag  coefficients 
of  the  Ith  mooring  part.  (CDN  and  CDT) 

Profile  Data; 

Depths,  greater  than  zero,  where  the 
current  changes  direction.  I goes  from 
1 to  IKK.  (Meters) 

For  SSD3.F0RT:  Relative  velocity  in 
the  (x,  y,  and  z)  direction  for  the 

Jth  zone.  J goes  from  1 to  IKK  + 1. 


(mm/sec) 
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For  SSD 3 1 . FORT : Absolute  velocity 
amplitude  in  the  1 ^ ‘ ' direction  for 
the  zone.  Following  data  are 

used  only  in  SSD31.FORT. 

WE:  Frequency  of  current  profile 

amp! i-udes  V(I,J).  (Rad/soc) 

O'”:  Time  step  for  calculation  of  the 

mooring  system  dynamic  equilibrium. 

T2:  Time  step  for  simulation  output. 

TMAX:  Maximum  tine  for  simulation. 

WS:  Surface  wave  frequency  (.).  (Rad/soc) 

AMP:  Surface  wave  amplitude  (.*  ) . (mm) 

BETA:  Surface  wave  direction  (;'■).  (Raaians) 

All  of  input  data  are  read  in  the  main  program. 
The  main  program  calls  subroutine  MOTION  which  calculates 
the  dynamic  equilibrium  of  the  subsurface  mooring  svsten  at 
discrete  times.  Subroutine  MOTION  calls  subroutine  FORCES, 
whenever  an  instrument  package  or  a subsurface  float  is 
encountered  in  the  mooring  line  incremental  integration 
scheme.  Both  MOTION  and  FORCES  call  subroutine  SPEED  to 
find  the  current  for  viscous  drag  calculations. 


5.0  CASE  STUDIES/SIMULATIONS 

This  section  presents  computer  simulations  of 
specific  case  studies;  using  the  computer  programs  outlined 
in  Section  4.  All  simulations  are  done  in  the  time-domain. 
Simulations  of  the  freely  floating  spar  buoy,  as  presented 
in  Section  5.1,  were  combined  to  obtain  the  response 
amplitude  operator  (RAO)  vs  Kh  (wave  number  multiplied  bv 
the  appropriate  buoy  draft).  Section  5.2  presents  the 
subsurface  moored  system  simulations.  Section  5.3  gives 
simulations  of  a surface  moored  system,  and  Section  5.4 
presents  simulations  of  a drifting  drogued  buoy. 

5.1  Spar  Buoy 

Both  the  two-dimensional  computer  program 
(SD2.FORT)  and  the  three-dimensional  computer  program 
(SD3.FORT)  are  used  to  simulate  each  of  the  two  spar  buoys; 
one  cylindrical,  and  the  other  tuned  presented  in 
Sections  5.1.1  and  5.1.2  respectively.  RAO  vs  Kh  plots 
a^e  presented  for  the  two-dimensional  simulations  of 

cylindrical  and  tuned  spar  buoys. 

i 

5.1.1  Cylindrical  Spar 

The  spar  buoy  (Figure  5.1)  used  in  these 

\ 

simulations  was  tested  in  a series  of  wave  tank  tests  to 
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derive  empirical  response  data  for  comparison  with  the 
computer  simulations.  This  comparison  has  not  been 
presented  in  this  report.  The  dimensions 
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Total  Weight  of  Buoy  = 

157.6  grams  (0.35  pounds) 

Buoy  Radius  of  Gyration  = 


Figure  5.1  Scale  Model  Spar  Buoy 
Description 


and  characteristics  of  this  scale  model  spar  buoy  are 
shown  in  Figure  5.1. 

Two  Dimensional  Simulations 

Figures  5.2,  5.3,  and  5.4  present  the  typical  surge, 
heave,  and  pitch  response  of  the  buoy  to  surface  waves  of 


rj'S'SS'S 
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Figure  5.2  Surge  Motions  of  the  Freely  floating  Cylindrical 
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amplitude  0.05  ft.  and  different  frequencies.  Table  5.1  is 
the  input  data  used  in  conjunction  with  the  computer  program 
SD2.FORT  to  obtain  these  simulations.  Noteworthy  is  the  fact 
that  the  viscous  (CDL)  and  pressure  (CDP)  drag  coefficients 
in  this  study  equal  zero.  Surge  drift  evident  in  Fig.  5.2 
is  due  to  the  bias  or  rectification  of  wave  exciting  forces, 
which  would  be  present  for  example  if  an  oscillatory  hydro- 
dynamic  force  were  to  act  on  a heaving  body  whose  motions 
were  at  varying  with  this  force.  Simulations  for  the  same 
data  as  Table  5.1  and  viscous  and  pressure  drag  coefficients 
equal  to  1.0  are  shown  in  Figures  5.5,  5.6,  and  5.7. 
Comparison  between  the  surge  motions  of  Figures  5.2  and  5.5 
shows  the  effect  of  viscous  drag  or.  surge  drift.  Surge  drift 
due  to  viscous  drag  occurs  as  the  submerged  area  of  the 
spar  and  its  inclination  vary  over  a wave  period,  and  is  a 
function  of  their  phase  relationships.  This  drift  is  most 
predominant  near  resonance.  In  this  case  heave  resonance 
occurs  at  Kh  = 1.0,  and  pitch  resonance  occurs  at  Kh  = 0.4. 
Heave  motions  (Figures  5.3  and  5.6)  are  not  effected  much 
by  viscous  drag  because  the  tangential  drag  due  to  CDL 
and  CDP  is  comparitively  small.  Comparison  of  Figures  5.4 
and  5.7  shows  that  the  system  frequency  present  in  the  top 
four  curves  of  Figure  5.4,  due  to  approximate  initial 
conditions,  is  damped  out  in  Figure  5.7  due  to  viscous 
drag  damping.  Otherwise,  the  amplitudes  do  not  show 
much  change. 
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Table  5.1  Input  Data  for  the  2-D  Analysis  of  the  Cylindrical  Spar  Buoy 


Figure  5.5  Surge  Motions  of  the  Freely  Floating  Cylindrical  Spar  Buoy 


Figure  5.6  Heave  Motions  of  the  Freely  Floating  Cylindrical  Spar  Bi  oy 
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Response  amplitudes  of  surge,  heave,  and  pitch; 

and  surge  drift  as  calculated  from  the  steady  state  portions 

of  the  times  histories  presented  in  Figures  5.2  through 

5.7,  and  other  similar  simulations  are  tabulated  in 

Table  5.2.  In  this  table  column  1 tabulates  the  frequencies 

(to)  of  waves  studied.  From  these  frequencies  Kh  is 

calculated  and  tabulated  in  column  2.  Columns  3 through  6 

are  four  columns  of  surge  amplitude  divided  by  £Q.  The 

first  of  these  four  columns  presents  the  results  of  the 

undamped  linear  analytical  solution  as  derived  in  Section  3.2. 

The  second  piesents  the  computer  simulation  results  with  £ 

equal  to  0.05  ft.  and  no  viscous  drag.  Third  and  fourth 

present  the  computer  simulation  results  with  viscous  drag 

and  £o  equal  to  0.05  ft.  and  0.25  ft.  respectively.  The 

simulated  plots  are  shown  for  two  wave  amplitudes  to 

emphasize  the  nonlinear  response  obtained  by  the  nonlinear 

theory.  Similarly  columns  7 through  10  tabulate  four  columns 

of  heave  amplitude  divided  by  £Q,  and  columns  11  through  14 

tabulate  four  columns  of  pitch  amplitude  divided  by  K£Q.  KCQ 

is  the  maximum  wave  slope.  Analytical  solution  for  surge 

drift  was  not  found  hence  only  three  columns  15  through  17 

2 2 

are  given  for  surge  drift  divided  by  Ko)£o  . Kw^0  is  the 
mean  "stokes  drift"  velocity  at  the  mean  free  surface.  Data 
points  of  Table  5.2  are  plotted  in  Figures  5.C  through  5.11 
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Analytical  Solution  - 1 


Figure  5.8  Surge  Magnification  Vs  Kh 
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Figure  5.9  Surge  Drift  Vs  Kh 
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to  present  plots  of  RAO  vs  Kh.  Figure  5.8  shows  the  buoy 

surge  per  foot  of  wave  amplitude  as  a function  of  Kh.  For 

zero  frequency  this  ratio  is  one  and  for  increasing 

frequencies  it  decreases  to  zero.  Figure  5.9  shows  the 

2 

buoy  drift  (in  surge  direction)  per  KwCQ  as  a function  of 
Kh.  Figure  5.10  shows  the  buoy  heave  per  foot  of  wave 
amplitude  as  a function  of  Kh.  This  ratio  is  one  for  zero 
frequency  and  exhibits  resonance  at  Kh  of  approximately  1.0. 
The  ratio  decreases  to  zero  with  increasing  frequencies. 
Figure  5.11  shows  the  buoy  pitch  per  K£o  as  a function  of  Kh. 
This  ratio  is  one  at  zero  frequency  and,  after  exhibiting 
resonance  phenomenon  at  around  0.4  Kh,  the  ratio  goes  to 
zero  with  increasing  frequencies.  Figures  5.8  through 
5.11  are  further  discussed  in  Section  6. 

Three  Dimensional  Simulations 

Figures  5.12  through  5.17  present  the  surge,  sway, 
heave,  roll,  pitch  and  yaw  responses  of  the  buoy  to 
surface  waves  of  the  same  amplitude  and  frequencies  as  in 
Figures  5.5  and  5.7.  Table  5.3  is  the  input  data  used  in 
conjunction  with  the  computer  program  SD3.FORT  to  obtain 
these  simulations.  Here  the  angle  8,  giving  the  wave 
propagation  direction  is  zero.  These  simulations  should 
give  the  same  response  as  Figures  5.5  through  5.7.  Comparing 
Figures  5.5  and  5.12,  we  see  a lot  of  difference  in  the 
surge  drifts.  Figure  5.12  looks  more  like  Figure  5.2  which 
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was  obtained  by  neglecting  viscous  drag.  Hence  we  see 

that  the  viscous  drag  effects  on  surge  drift  which  show 

up  in  the  two-dimensional  non-linear  simulations  (Figure  5.5), 

do  not  show  up  in  the  three-dimensional  simulation 

(Figure  5.12)  where  assumptions  of  small  motions  are  made. 

This  assumption  may  be  in  error  the  most  near  the 
resonance  frequency  (Heave  resonance;  Kh  = 1.0.  Pitch 
resonance;  Kh  = 0.4),  where  the  two  Figures5.5  and  5.12 
differ  the  most.  Comparison  of  Figures  5.14  and  5.16 

to  Figures  5.6  and  5.7  shows  only  slight  differences. 

Figures  5.13,  5.15,  and  5.17  show  no  response  as  expected. 
Simulations  for  the  same  data  as  Table  5.3  but  g = n/4 
are  presented  in  Figures  5.18  through  5.23.  These  figures 
show  the  resolution  of  surge  in  Figure  5.12  to  surge  and 
sway  in  Figures  5.18  and  5.19;  and  the  resolution  of  pitch 
in  Figure  5.16  to  roll  and  pitch  in  Figures  5.21  and  5.22. 

Yaw  response  remains  the  same  as  expected. 

5.1.2  Tuned  Spar 

The  spar  buoy  used  in  these  simulations  along 
with  its  dimensions  and  characteristics  is  shown  in 
Figure  5.24. 

Two-Dimensional  Simulations 

Figures  5.25,  5.26,  and  5.27  show  the  typical  surge, 
heave,  and  pitch  response  of  the  buoy  to  surface  waves  of 
amplitude  0.5  ft.  and  different  frequencies.  Table  5.4 
is  the  input  data  used  in  conjunction  with  the  computer 
program  SD2.FORT  to  obtain  these  simulations.  Here  again 
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viscous  and  pressure  drag  coefficients  equal  zero. 

Figure  5.28  through  5.30  use  the  same  data  as  Table  5.4 
except  that  viscous  and  pressure  drag  coefficients  are 
equal  to  1.0.  Once  again  comparisons  between  the  two  sets 
of  figures  (5.25  through  5.30)  show  similar  differences 
as  explained  for  the  cylindrical  spar  buoy.  In  this  case 
the  heave  resonance  occurs  at  Kh  = 0.22  and  the  pitch 
resonance  at  Kh  = 0.4.  Also  in  this  case  the  viscous 
damping  in  heave  and  pitch  responses  is  more  pronounced. 

Response  amplitudes  of  surge,  heave,  and  pitch; 
and  surge  drift  as  calculated  from  the  steady-state  portions 
of  the  time  histories  presented  in  Figures  5.25  through  5.30, 

and  other  similar  simulations  are  tabulated  in  Table  5.5 
which  is  similar  to  Table  5.2.  In  this  case  the  wave 

j 

amplitudes  studied  are  0.5  ft.  and  2.5  fc.  respectively. 

Data  points  of  Table  5.5  are  plotted  in  Figures  5.31 

through  5.34  to  presents  plots  of  RAO  vs  Kh.  These 

plots  are  similar  to  Figures  5.8  through  5.11.  Figure  5.33b 

is  an  expanded  view  of  Figure  5.33a,  near  the  heave 
resonance  frequency.  The  ratio  between  buoy  heave  and  the 

wave  amplitude  is  one  at  zero  frequency  and  exhibits  resoance 
at  Kh  of  approximately  0.2.  The  ratio  decreases  to  zero  at  Kh 
equal  to  0.32,  after  which  it  increases  a little  before 
decreasing  to  zero  with  increasing  frequencies.  Such  behavior, 
where  the  neave  response  is  minimized  at  a given  freouency  is 
the  main  function  of  a tuned  (as  opposed  to  cylindrical)  spar 
buoy. 
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Table  5.5  Dynamic  Response  (Surge,  Heave  and  Pitch  Amplitude; 
and  Surge  Drift)  of  a Tuned  Spar  Buoy  to  Waves 


Figure  5.32  Surge  Drift  Vs  Kh 
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Three-Dimensional  Simulations 
Figures  5.35  through  5.40  present  the  surge, 
sway,  heave,  roll,  pitch,  and  yaw  responses  of  the  buoy 
to  surface  waves  of  the  same  amplitude  and  frequencies  as 
in  Figures  5.28  to  5.30.  Table  5.6  is  the  input  data 
used  in  conjunction  with  the  computer  program  SD3.F0RT 
to  obtain  these  simulations.  Here  8 equals  zero. 

Differences  between  Figure  5.28  and  5.35  can  be  attributed 
directly  to  the  assumptions  of  three-dimensional 
analysis.  Simulations  for  the  same  data  as  Table  5.6 
and  8 = fl/4  are  presented  in  Figures  5.41  through  5.46. 

5 . 2 Subsurface  Moored  System 

Three  dimensional  computer  programs  (S3S3.FORT, 
SSD3.FORT  and  SSD31.FORT)  described  in  Section  4.2  are 
used  for  these  simulations.  Mooring  system  used  as  a 
case  study  for  these  simulations  is  shown  in  Figure  5.47. 

This  mooring  system  was  actually  deployed  during  the  MODE 
experiment  and  was  called  mooring  No.  1 (station  481) . 

Figures  5.48  and  5.49  display  the  three  coordinates 
(x,  y,  z)  of  the  top  of  this  mooring  line  in  response  to 
the  relative  velocity  data  displayed  in  Figures  5.50  and 
5.51.  This  simulation  is  a part  of  the  study  (CHHABRA,  1976), 
where  the  current  record  from  the  topmost  vector-averaging 
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Figure  5.49  trajectory  of  the  Subsurface  I-'oored  System  'Ibp 
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current  meter  of  this  mooring  was  corrected  for  the  effects 
of  mooring  motion,  and  power  spectra  of  the  uncorrected 
and  corrected  signals  were  compared.  The  relative  velocity 
data  as  shown  in  Figures  5.50  and  5.51  was  sampled  every 
fifteen  minutes  at  eight  locations  on  this  moori'  line. 
Computer  program  SSS3.FORT  was  used  for  this  simulation 
and  no  surface  wave  was  assumed  present.  Figures  5.52 
and  5.53  show  this  same  response  with  the  computer  program 
SSD3.FORT,  which  calculates  and  includes  inertia  forces. 
Overlay  of  the  two  responses  shows  no  noticeable  difference. 
Input  data  used  with  these  two  computer  programs,  describing 
the  system  are  shown  in  Tables  5.7  and  5.8. 

Response  of  the  top  of  this  mooring  line  to 
hypothetical  absolute  velocity  data  (Figures  5.54  and  5.55), 
is  shown  in  Figures  5.56  and  5.57.  This  simulation  was 
done  with  computer  program  SSD31 . FORT  and  the  input  data 
shown  in  Table  5.9.  Surface  wave  was  again  assumed  to  be 
absent  in  this  simulation. 

5. 3 Surface  t.oored  System 

The  three-dimensional  computer  program  SD3.IORT 
is  used  for  this  simulation.  The  surface  moored  system 
simulated  is  shown  in  Figure  5.58.  This  system  of  spar 
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Table  5.7  (cent.)  Input  Data  for  the  Subsurface  floored  System 
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buoy/buoyant  tether/subsurface  mooring  line  with  an 
instrument  line  hanging  beneath  the  spar  was  deployed 
as  part  of  the  ONR/NDBO  Mooring  Dynamics  Experiment 
during  October  1976.  The  system  is  reduced  to  seven 
nodes  for  this  simulation.  These  nodes  are  shown  in 
Figure  5.59.  Figure  5.59  shows  the  mean  configuration 
of  these  nodes  in  the  x-z  plane.  This  configuration  is 
the  result  of  the  static  solution  calculated  by  SD3.F0RT 
using  the  equivalent  velocity  profile  described  in 
Section  4.1.  The  equivalent  velocity  profile  and  the 
actual  velocity  profile  in  the  x-z  plane  is  also  shown  in 
Figure  5.59.  Figures  5.60  through  5.64  display  time  responses 
of  these  seven  nodes  to  the  input  data  presented  in  Table  5.10 
which  describes  the' moored  system  and  the  environment  forcing 
this  system.  Forcing  the  system  is  a constant  current 
profile  and  a surface  wave  of  amplitude  1 foot,  a frequency 
of  1 radian  per  second,  and  having  a wave  direction 
inclined  at  0.1  radians  with  the  x-axis  of  the  earth 
fixed  frame.  Figures  5.60  through  5.62  are  the  plots  of 
surge,  sway  and  heave  responses  of  the  seven  nodes. 

Figure  5.63  displays  the  roll,  pitch  and  yaw  of  the  tuned 
spar  buoy  (node  6),  and  Figure  5.64  displays  tensions  in 
the  seven  segments  of  mooring  lines  preceding  the  seven  nodes. 


Figure  5.61  Sway  J lotions  of  a Surface  Moored  Systran 
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5 . 4 Drifting  Droguol  Duoy  System 

A submerged  window  shade  drogue  attached  to  a 

cylindrical  spar  buoy  by  means  of  an  elastic  nylon  line 

is  simulated  on  the  computer  by  the  three  dimensional 

computer  program  SD3.F0RT.  The  system  is  subjected  to 

simple  harmonic  surface  waves  without  any  current  profile. 

A 4"?  feet  long  3/S"  nylon  line  is  attached  to  the  center 

of  the  bottom  of  the  buoy.  The  spar  buoy  is  30  feet  long 

and  0.66  foot  in  diameter.  Its  static  draft  with  no  lines 

attached  is  19.43  feet  and  weighs  426  pounds  in  air. 

The  distance  between  the  center  of  gravity  and  the  bottom 

end  of  the  buoy  is  8.69  feet.  The  window  shade  drogue 

2 

is  made  of  nylon  cloth  (9  1/4  oz/yd  ) and  has  the 
dimensions  of  7.5'  x 32'.  Its  approximate  thickness  is 
0.02".  The  drogue  has  two  metallic  bars,  one  at  each  end. 

The  bottom  bar  also  acts  as  the  dead  weight,  keeping  the 
drogue  in  tension.  For  computer  simulation  the  drogue 
was  divided  into  four  strips  and  the  parameters  of  each 
strip  lumped  at  a node.  Tensions  in  three  links  between 
these  four  nodes  and  (he  tension  in  the  nylon  tether  line 
were  simulated  on  the  computer.  Figure  5.6r>  shows  the 
system.  Figures  5.66  through  5.70  show  the  tine  histories 
of  the  response  of  the  system  in  a surface  wave  of 
amplitude  1 foot  and  circular  frequency  of  1 radian  per 
second,  and  having  a wave  direction  inclined  at  0.25  radians 
with  the  x-axis.  No  current  profile  is  present.  Figures  5.66 
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Figure  5.65  Drifting  Drogued  Buoy  System 
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through  5.6  8 show  the  surge,  sway  and  heave  coordinates 
of  the  spar  bu<.y  C.G.  and  the  four  nodes  on  the  drogue. 
Figure  5.69  shows  the  roll,  pitch,  and  yaw  angles  of  the 
buoy  and  Figure  5.70  shows  the  tension  magnitudes  in 
the  nylon  line  helow  the  buoy  and  the  three  segments  cf 
the  drogue.  The  input  data  used  for  this  simulation  is 
presented  in  Table  5.11. 
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6 . 0 SUMMARY 

Dynamics  of  moored  and  drifting  buoy  systems  in 
a three-dimensional  space  has  been  presented  in  this 
report.  Formulation  of  mathematical  models,  which  are 
programmed  on  the  inhouse  computer  AMDAHL  <70  V6,  along 
with  time-domain  computer  simulations  of  these  systems  are 
presented.  Four  case  studies  {freely  floating  spar  buoy; 
a single  point  subsurface  moored  system;  a spar  buoy  plus 
instrument  line/buoyant  tether/subsurface  mooring;  a 
freely  drifting  spar  buoy  attached  to  a window  shade  drogue) 
are  simulated  with  forcing  being  supplied  by  a velocity 
profile  and  a fully  developed  surface  wave  field.  One 
(a  single  point  subsurface  moored  system)  of  these  four 
case  studies  was  evaluated  with  full  scale  ocean  test 
data  (CHHABRA,  DAHLFN  and  FROIDEVAUX,  1974;  CHHAERA,  1976). 
Simulations  of  other  three  case  studies  can  now  be  readily 
evaluated  with  full  scale  ocean  *rest  data. 

For  the  two-dimensional  analysis  of  a freely 
floating  spar  buoy,  where  the  assumption  of  small  motions 
is  not  made,  three  sets  of  computer  simulations  are 
compared  with  the  undamped  linear  analytical  solution. 

, 

IThis  comparison  is  presented  in  Figures  5.3  through  5.11 
for  a cylindrical  spar  and  Figures  5.31  through  5.34 
for  a tuned  spar  buoy.  This  comparison  shows  the  effects 
of  viscous  drag  and  varying  wave  amplitude  on  the  response 
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of  a freely  floating  spar  buoy.  Plots  for  the  undamped 
linear  analytical  solution  and  the  computer  simulation 
where  wave  amplitude  is  small  and  viscous  drag  is  neglected 
are  almosT.  identical.  For  the  case  of  a cylindrical  spar, 
when  viscous  drag  is  added,  the  surge  magnification  vs 
Kh  plots  show  a scatter  near  the  heave  and  pitch  natural 
frequencies  (Figure  5.8).  This  emphasizes  the  coupling 
effects.  Surge  drift  vs.  Kh  plots  in  Figure  5.9  show 
a markedly  different  responses  near  these  same  frequencies 
and  viscous  damping  is  evident  in  heave  and  pitch 
magnification  plots  (Figures  5.10  and  5.11).  The  non- 
linearity in  the  response  with  respect  to  wave  amplitude 
is  also  evident.  These  plots  also  show  a heave-pitch 
coupling  near  their  resonance  frequencies.  Similar 
observations  can  be  made  for  the  case  of  the  tuned  spar. 

For  the  three-dimensional  analysis  of  a freely 
floating  spar  buoy,  comparison  is  made  with  the  two- 
dimensional  analysis  for  the  same  input  data.  The  only 
appreciable  difference  is  observed  ir.  the  surge  drift, 
which  can  be  attributed  directly  to  the  small  motion  assumption 
of  the  three-dimensional  analysis. 

Simulation  of  the  spar  buoy  plus  instrument 
line  tethered  to  a subsurface  mooring  line  shows  a note- 
worthy transmittal  of  motion  from  the  surface  buoy  to  the 
subsurface  line,  by  the  tether.  The  tether  line  seems  to 
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move  mainly  in  the  tangential  direction  due  to  the  high 
viscous  drag  force  in  the  transverse  direction.  As  a 
result  the  combined  surge,  sway,  and  heave  motions  of  the 
spar  are  transmitted  to  the  subsurface  line.  Hence,  the 
surge/sway  motions  of  the  spar  (even  if  heave  was 
negligible)  could  be  transmitted  to  the  subsurface  line  as 
heave.  Therefore  some  modifications  of  the  system  may  be 
necessary  to  minimize  longitudinal  motions  of  the  subsurface 
line . 

Simulation  of  the  freely  drifting  spar  buoy 
attached  to  a window  shade  drogue  shows  drift  of  the  buoy 
and  the  drogue  even  though  no  current  (except  surface  wave) 
is  present-  This  drift  is  due  to  the  same  reasons  as 
postulated  for  the  freely  floating  spar  buoy. 

Full  scale  ocean  test  data  for  the  spar  buoy  plus 
instrument  line/buoyant  tether/subsurface  mooring  is  now 
available  (October  1976  0N3/NDB0  Mooring  Dynamics  Experiment) . 
Using  this  and  other  data  which  might  be  available  in  the 
future;  the  next  logical  step  is  the  evaluation  of  the 
mathematical  models  of  the  remaining  three  case  studies. 
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APPENDIX  A 

Computer  Program  Listings  Pago 

SD3.F0RT  197 

SD2.F0RT  222 

SSD31 . FORT  245 

SSS  3 . FORT  255 


SSD3.F0RT 


262 


SD3.F0RT 
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